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Contemporary understanding of correlations in quantum many-body systems and in quantum phase transi¬ 
tions is based to a large extent on the recent intensive studies of entanglement in many-body systems. In con¬ 
trast, much less is known about the role of quantum nonlocality in these systems, mostly because the available 
multipartite Bell inequalities involve high-order correlations among many particles, which are hard to access 
theoretically, and even harder experimentally. Standard, ’’theorist- and experimentalist-friendly” many-body 
observables involve correlations among only few (one, two, rarely three...) particles. Typically, there is no 
multipartite Bell inequality for this scenario based on such low-order correlations. Recently, however, we have 
succeeded in constructing multipartite Bell inequalities that involve two- and one-body correlations only, and 
showed how they revealed the nonlocality in many-body systems relevant for nuclear and atomic physics [Sci¬ 
ence 344, 1256 (2014)]. With the present contribution we continue our work on this problem. On the one hand, 
we present a detailed derivation of the above Bell inequalities, pertaining to permutation symmetry among the 
involved parties. On the other hand, we present a couple of new results concerning such Bell inequalities. First, 
we characterize their tightness. We then discuss maximal quantum violations of these inequalities in the general 
case, and their scaling with the number of parties. Moreover, we provide new classes of two-body Bell in¬ 
equalities which reveal nonlocality of the Dicke states—ground states of physically relevant and experimentally 
realizable Hamiltonians. Finally, we shortly discuss various scenarios for nonlocality detection in mesoscopic 
systems of trapped ions or atoms, and by atoms trapped in the vicinity of designed nanostructures. 


I. INTRODUCTION 


Nonlocality is a property of correlations that go beyond the paradigm of local realism [1-5]. According to the celebrated 
theorem of J. S. Bell, correlations between the results of measurements of local observables on some quantum states might be 
nonlocal, i.e. they do not admit a local hidden variables (LHV) model a la Einstein-Podolski-Rosen (EPR) [1,6], In such 
instances we talk about Bell nonlocality. Operationally, such correlations cannot be simulated by correlated classical random 
variables. This can be detected by means of the so-called Bell inequalities [2] - the celebrated example of such is the famous 
Clauser-Horne-Shimony-Holt inequality (CHSH) [7]. Bell nonlocality is interesting for at least three reasons: 

• It is a resource for certain information-processing tasks, such as secure key distribution [8-10] or certified quantum 
randomness generation [11-13]. Hence, revealing the nonlocality of a given composite quantum state is one of the central 
problems of quantum information theory, of relevance for future quantum technologies. 

• It lies at the heart of fundamental aspects of quantum physics [4, 6, 14], e.g., playing a fundamental role in our under¬ 
standing of intrinsic quantum randomness, and sometimes even in philosophical contexts such as the free will problem 
[15, 16], 

• It is extremely challenging mathematically: characterization of non-locality is a complex and difficult problem. In fact it 
was proven that determining whether a given set of correlations belongs or not to the set of local EPR correlations, is NP- 
complete, whereas specifying all Bell inequalities for growing number of involved parties, and/or measured observables, 
and/or outcomes of measurements is NP-hard ([17, 18]; see also [3, 5] and references therein). 

Quantum states that violate a Bell inequality are necessarily entangled, or in other words, they are not separable and cannot 
be represented as mixtures of projections onto product states [19] (for a review on entanglement see [20]); the opposite is not 
necessarily true. Already in 1991 Gisin proved [22] that any entangled pure state of two parties violates a Bell inequality (in 
fact, they violate the CHSH inequality). This result was extended to pure entangled states of an arbitrary number of parties by 
Popescu and Rohrlich [23]. However, the situation changes completely when moving to mixed states: In his seminal paper from 
1989 [19] Werner constructed examples of mixed bipartite states that admit a LHV model for all local projective measurements, 
and nevertheless are entangled. This result was then generalized by Barrett to arbitrary measurements [26], Very recently, some 
of us have shown that entanglement and nonlocality are inequivalent for any number of parties, i.e. that there exist genuinely 
multipartite entangled states that do not display genuinely multipartite nonlocality [27]. It is then clear that entanglement is a 
weaker property of quantum states than nonlocality. 

Over the last decade, entanglement has proven to be very useful to characterize properties of many-body systems and the 
nature of quantum phase transitions (QPT) [28]. Let us illustrate this statement for lattice spins models described by local 
Hamiltonians, that is, those involving finite range interactions, or interactions rapidly decaying with distance. Most of the results 
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listed below also hold for lattice Bose and Fermi models, and even for quantum field theories. In particular, in the ground states 
(GS) (low energy states) of such models the following properties are true (for general reviews see Refs. [29, 30]): 

• If one calculates the reduced density matrix for two spins, this density matrix typically exhibits entanglement for short 
separations of the spins only. This occurs even at criticality; still entanglement measures show signatures of QPTs [31, 32]; 

• If, instead of calculating a reduced density matrix, one tries to perform optimized measurements on the rest of the system, 
one concentrates then entanglement in the chosen two spins. One obtains in this way localizable entanglement [33]; the 
corresponding entanglement length diverges when the standard correlation length diverges, i.e., at standard QPTs. Inter¬ 
estingly, there exist critical systems for which the correlation length remains finite, while the entanglement localization 
length goes to infinity [34]; 

• For non-critical systems, GSs and low energy states exhibit the, so-called, area laws: the von Neuman (or Renyi) entropy 
of the reduced density matrix of a block of size I scales as the size of the boundary of the block, dl\ at criticality a 
logarithmic divergence occurs frequently [35] (for a review see [36, 37]). These results are very well established in ID, 
while there are plenty of open questions in 2D and higher dimensions; 

• GSs and low energy states away from criticality can be efficiently described by the, so-called, matrix product states, or 
more generally tensor network states (cf. [38]); even at criticality MPS are efficient in ID. 

• Topological order (at least for gapped systems in ID and 2D) exhibits itself frequently in the properties of the, so-called, 
entanglement spectrum, i.e., the spectrum of the logarithm of the reduced density matrix of a block I [39], and in 2D in 
the appearance of the, so-called, topological entropy, i.e., a negative constant correction to the area laws [40, 41]. 

A natural question thus arises: Does nonlocality also play an important role in characterizing correlations in many-body 
systems? Apart from its fundamental interest, so far the role of nonlocality in such systems has hardly been explored (see e.g. 
[21]). Entanglement and nonlocality are known to be inequivalent quantum resources. In principle, a generic many-body state, 
say a ground state of a local Hamiltonian, is pure, entangled and, because all pure entangled states violate a Bell inequality [23], 
it is also nonlocal. However, this result cannot be easily verified in experiments, because the known Bell inequalities (see, e.g., 
[42-46]) usually involve products of observables of all observers (also referred to as parties). Unfortunately, measurements of 
such observables, although in principle possible [47, 48], are technically extremely difficult; instead one has typically ’’easy” 
access to few-body correlations, say one- and two-body, in generic many-body systems. Thus, the physically relevant question 
concerning the nonlocality of many-body quantum states is whether its detection is possible using only two-body correlations. 

When tackling this problem one has to face the following technical challenges. First of all, finding Bell inequalities that 
are valid for an arbitrary number of parties is usually a difficult task, as the complexity of characterizing the set of classical 
correlations scales doubly exponentially with the number of parties. Then, it is natural to expect that Bell inequalities based on 
only few-body correlations are, in general, weaker in revealing nonlocality than those constructed from high-order correlations 
(that is, correlations among many of the observers). Whereas Bell inequalities involving correlations among all but one parties 
have been constructed in [49-51], thus proving that all-partite correlations are not necessary to reveal nonlocality, here we have to 
consider the more demanding question of whether nonlocality detection is possible for systems of an arbitrary number of parties 
from the minimal information achievable in a Bell test, i.e., two-body correlations. We answered this question positively in Ref. 
[52] by proposing classes of Bell inequalities constructed from one- and two-body expectation values, and, more importantly, 
providing examples of physically relevant many-body quantum states (i.e., GSs or low energy states of physically relevant 
Hamiltonians) that violate these inequalities. Note that finding and classifying such states is an interesting task in itself, especially 
in view of the fact that many genuinely entangled quantum many-body states have two-body reduced density matrices (or, in 
other words, covariance matrices) that are compatible with two-body reduced density matrices of some separable state. This is 
the case of the so-called graph states, as demonstrated in Ref. [53]; obviously one cannot detect entanglement of such states 
with two-body correlators, not even mentioning nonlocality. 

It should also be stressed that an analogous question was already explored in the case of entanglement, which next to nonlo¬ 
cality is a key resource quantum information theory. Several entanglement criteria relying solely on two-body expectation values 
have been proposed [54]. In particular, in [55] the possibility of addressing two-body statistics (although not individually) via 
collective observables was exploited. 

The main aims of this work are two-fold. First, we present a detailed derivation of multipartite Bell inequalities involving 
only one- and two-body correlators, and pertaining to permutation symmetry between the involved parties presented in [52]. 
To this end we characterize the polytope of LHV correlations in such case. We then discuss the maximal quantum violations 
of these inequalities in the general case, including its scaling with the number of parties. We pay particular attention to Bell 
inequalities violated by the family of Dicke states that are the ground states of physically relevant and experimentally realizable 
Hamiltonians. We discuss shortly various scenarios for detection of non-locality in mesoscopic systems of trapped ions or atoms, 
and by atoms trapped in the vicinity of designed nanostructures. The second aim is to strongly develop and generalize the results 
of Ref. [52], First, we impose conditions under which our Bell inequalities are tight. Then, we provide a new class of two-body 
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Bell inequalities which reveal nonlocality in all the entangled Dicke states, not only those studied in [52], Finally, we construct 
and characterize an analytical class of pure states that well approximate states maximally violating some of our Bell inequalities. 
With their aid we can analytically compute the maximal violation of these Bell inequalities in the thermodynamic limit. 

The manuscript is organized in sections, in which we define the problem and present the results. All the technical proofs and 
demonstrations, as well as a discussion of concrete experiments aimed at generation of the Dicke states, are shifted to the three 
Appendices. In Section II we describe the generic Bell experiment, using the contemporary language of device-independent 
quantum information theory in terms of conditional probabilities for n parties compatible with the no-signalling principle. The 
structure of the polytope defined by the local correlations (local Bell polytope) is discussed in Section III. Here we first describe 
the general case, then how we can restrict the analysis to a smaller polytope using a general symmetry group, and finally we 
focus on the permutationally invariant Bell polytope. In Section IV various classes of Bell inequalities are derived and discussed. 
Section V is devoted to results concerning maximal quantum violation of the derived inequalities. We present here also analytical 
expressions for a class of states that violate the derived Bell inequalities, and such that the analytical violation converges to the 
maximal numerical one as n becomes large. Various aspects of the robustness of these inequalities and their violations are 
analyzed here, with an emphasis on the limit of a large number of parties. In Section VI we focus on inequalities detecting 
nonlocality of the entangled Dicke states, and discuss also their parameter dependence. Section VII contains a short conclusion, 
and a discussion of several promising experimental systems in which the proposed Bell test could be successfully implemented: 
mesoscopic systems of ultracold trapped ions, neutral atoms in a single, or double well trap, or atoms near tapered optical fibers 
and nanostructured optical crystals. 

Appendix A contains the proofs of all the theorems from Sections III-VI, while Appendix B presents some details concerning 
the block-diagonalization of the Bell operator, necessary to determine the maximal quantum violation. As mentioned above. Ap¬ 
pendix C contains detailed calculations concerning the robustness to experimental errors of the 2-body Bell inequalities. These 
calculations are relevant for the potential implementation of the derived Bell inequalities, for instance in the recent experiment 
of Ref. [56], 


II. THE BELL EXPERIMENT 

Let us consider the standard multipartite Bell-type scenario in which n spatially separated observers, denoted Ai- ■ ■ ■, A„,_j 
share some //-partite resource. On their share of this resource every observer A, can perform one of to measurements, each 
having d outcomes. In what follows these measurements are denoted by A4x- with x, £ {0,..., m — 1}, while their outcomes 
by di £ {0,..., d — 1}. This scenario is usually referred to as (n, to, d). 

The correlations between results obtained in the above experiment are described in terms of the conditional probabilities 

P(a 0 ,...,a n _ i|x 0 ,... ,x n _i), 0 < a, < d, 0 < Xi < to, (1) 

(i) 

that the ith party obtained a, upon performing the measurement M. Xi - 

It is convenient to collect all the probabilities in (1) into a vector with ( md) n components, each corresponding to a different 
combination of outcomes a o, ■ ■ •, i and measurements choices xo,..., x ra -i- Let us call such vector p. Clearly, one must 

ask for probabilities to be non-negative, that is, P(ao,..., a n -i\xQ ,..., x„_i) > 0 must hold for all choices of ao,..., a n _i 
and too, - - -, x n -i, and normalized, meaning that 

^ ^ P(o 0 , ■ ■ ■; Oj n — 11tco? * • • • tr n _r) — 1 (2) 

a 0v-5 a n-l 

is satisfied for any sequence xo, ■ ■ ■, x„_i ■ The latter condition implies that only m n (d n — 1) of the components of p are 
independent. Let us also notice that all the above constraints define a region in space, which is a polytope, denoted Pg, to which 
any mathematically consistent p belongs. Recall that a polytope is a compact convex set with a finite number of extreme points. 

Now, depending on the nature of the n-partite resource that the parties share, different types of correlations can be obtained. 
We are mainly interested in three scenarios: the case in which p admits a local hidden variable (LHV) model, the case in which p 
can be obtained through performing local measurements on quantum states (Q) and the case in which p fulfils the no-signalling 
principle (NS). 

LHV correlations. Those conditional probability distributions that p have a local hidden variable model have the following 
form: 


P(&0? * • • ) tZ n _i |x 0 , • • • , X n — l) 


« n —1 

/ pw n P(a,i\xi, A)dA, 

4=0 


(3) 


where A £ A is some hidden variable, distributed according to a probability density p{ A), and A is the space of hidden variables. 
Operationally, the interpretation of (3) is that the outcome a; of the x,-th measurement is produced locally at site i as a function 
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of the input Xi and the value of this hidden variable A. The resulting probability distribution is the one corresponding to our 
understanding of the classical correlations between local independent measurements: In order to simulate correlations p of this 
form, parties could gather beforehand and agree on which A to pick at each run of the experiment, distributed according to 
p( A). This information, called also shared randomness is precisely the resource needed—given A and aeach party would just 
simulate P(a,i\xi, A) on their own. The feasibility region for LHV model probability distributions is a polytope and its vertices 
are given by those p in which all parties follow a local deterministic strategy. We shall denote this polytope P. 

The vertices of P are straightforward to construct as these are those p’s for which all probabilities (1) factorize: 

n— 1 

P(a 0 5 • • * ) On— 1 ko i - 5 ^n— l) 11 P(0 J i |a^) , (4) 

i=0 

and each factor is deterministic, i.e., P(ai\xi) = S(ai = a? 4 ), where S is the Rronecker delta function and 0 < a < d is a 
different integer for every vertex. One immediately finds that in the (n, to, d) scenario the number of vertices is d" vn . 

Quantum correlations. Imagine now that the resource the parties share is some n-partite quantum state p acting on a product 

Hilbert space T-L n = (C D )®" of local dimension D. In such case, each local measurement A4^'; can be thought of as a collection 

( x .) 

of positive semi-definite operators H K ai corresponding to outcomes a* = 0 ,..., d — 1 and obeying 

d -1 

E n i? } = ^ (5) 

CLi —0 

for every x r and i. Then, the probabilities (1) can be represented via Born’s rule: 

P(a 0 ,...,a n - i|x 0 ,...,x„_i) = TV [p (g ^ 1 ) . ( 6 ) 

The feasibility region for probability distributions of the form ( 6 ) is a convex set, denoted Q. Unlike the local set, Q is not a 
polytope and hence its boundary is more difficult to characterize. Still, it can be well approximated [57], 

Non-signalling correlations. For completeness of this work it is worth recalling the notion of non-signalling correlations. 
These are those p’s that fulfil the no-signalling principle which says that information cannot be transmitted instantaneously. This 
means that in the above Bell scenario the choice of measurement made by one party cannot influence the statistics seen by the 
remaining parties. In terms of the conditional probabilities (1) this is equivalent to a set of linear constraints of the form 

E P(a 0 , ...,a l ,... a n -i\xo, ...,x i} ..., x n -i) = E P ( a o> a„_ i |tr 0 ,..., x\,..., x n -{). (7) 

CLi CLi 

for all Xi ^ x' i and ao ,..., a^-i, aj+i,..., ad -i and Xo , ■ • •, Xi+i ,..., Xd-i and all i. The feasibility region correspond¬ 
ing to distributions p satisfying the no-signalling principle, which we denote P ns> is again a polytope, as it is an intersection 
of the polytope P 5 with the linear subspace defined by (7). However, unlike in the case of the local polytope P, the vertices of 
Pjvs are in general difficult to characterize, and except for the simplest scenarios ( 2 , 2 , d ) [60] and ( 2 , to, 2 ) [ 111 ] and they are 
unknown (see also Ref. [59]). 

When the no-signalling constraints and the normalization conditions are taken into consideration, the number of independent 
variables of any p S P ns is reduced to [m(d — 1) + l] n — 1 which is the dimension of the polytope P ns an d so is the dimension 
of the remaining two sets P and Q. Importantly, these three sets are different. In fact, one has the following chain of inclusions: 

PCQCP WS C p s . ( 8 ) 

First, P C Q because all correlations belonging to P can be realized with separable states. Yet, due to the celebrated Bell’s 
theorem P 7 ^ Q [2, 5], Second, Q C P,v 5 because all correlations given by Eq. ( 6 ) trivially satisfy the no-signalling principle. 
On the other hand, there exist correlations obeying the no-signalling principle which are not quantum, with the most paradigmatic 
example being the so-called Popescu-Rohrlich box [23], Hence, Q -/ P ns- Finally, it is easy to verify that not all elements of 
Ps satisfy the no-singalling principle, and therefore P vs Q Ps- 

Due to the structure of the local polytope P, the natural way of checking whether p £ Q \ P is to use Bell inequalities. These 
are linear inequalities formulated in terms of the probabilities resulting from the local measurements performed by the observers, 
that is, 


/ :=^2 T a , x p {a\x) > ~/3c (9) 

a,x 

where T a x e IF! are some coefficients and Be = - min^ g p I is the so-called classical bound of the Bell inequality. Violation 
of such inequalities signals nonlocality. 
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All polytopes admit a dual description: They can be completely characterized either by listing all their vertices (i.e., a point p 
belongs to a poly tope Ps if and only if it is a convex combination of its vertices) or as the intersection of a number of half-spaces. 
Such intersection can be taken to be minimal and in this case every of these half-spaces intersected with P 5 defines a facet of 
Pg. Nevertheless, going from one description to the other is not a simple task in general. The best known algorithm [58] has 
complexity 0(v L-°/ 2 J), where v is the number of vertices (facets) of the polytope and D is the dimension of the space in which 
the polytope lives. In the case of P, its vertices, although there are many, are straightforward to find, but not its facets, which 
give rise to the so-called tight Bell inequalities. Only few scenarios have been completely characterized, none of them with more 
than 3 parties [44, 112]. In the case of Pjvs the situation is the opposite: its facets are easy to construct, since they are the ( md) n 
positivity constraints on every probability (1), but its vertices, known as extreme non-signalling boxes, are hard to find. In this 
case only few scenarios have been solved [26, 59, 111], 

For further benefits let us eventually mention that in the case d = 2, which is the case we mostly consider in this work, it is 
convenient to re-label outcomes from 0,1 to ±1 (i.e., a 1 —> (— 1 )“), and re-express all probabilities ( 1 ) in terms of expectation 
values 

“n 

for 0 < ij < ... < ifc < n, j 1 £ {0,1} and 1 < k < n. Below, the number of parties appearing in such a correlators will also 
be referred to as order of the correlator; for instance, correlators (M^M^) are of order two. As a direct consequence of (4), 
the vertices of P are in this case given by 

= ( 11 ) 


where each local expectation value (Mx}) equals either — 1 or 1. 


III. PROJECTING THE LOCAL POLYTOPE 

In the present work we are interested in the study of non-locality in many-body quantum systems consisting of a large number 
of parties n. Finding all the facets of the corresponding local polytope P (that give rise to the optimal Bell inequalities) would 
completely characterize it and would be the first step to answer the above question. Unfortunately, even in the simplest case of 
to = d = 2 but for large n, P is an object of tremendous complexity, and attempting to characterize it following this approach 
becomes an intractable task. Since both its dimension and its number of vertices are exponential in n, we aim at projecting P 
onto a simpler object, which has less vertices and which is embedded in a lower-dimensional space. 

This simplification does not come for free and the obtained Bell inequalities would, in general, be weaker. If one applies 
a projection n to P, then 7r(P) can only give sufficient conditions for nonlocality [45]. If some correlations p are nonlocal, 
i.e., p (ji P, it may happen that it(p) ^ 7r(P) or 7 r(p) £ 7r(P), depending on the projection 7r. However, if 7 t(p) ^ 7r(P), 
then p ^ P. Hence, violating a Bell inequality corresponding to 7r(P) certifies nonlocality, while proving that correlations 
that lie inside 7r(P) is inconclusive. In this section we then discuss how to choose 7r in such a way that, on one hand, the 
Bell inequalities constraining 7r(P) are powerful enough to reveal nonlocality for arbitrary n, and, on the other hand, they are 
accessible with current experimental technology. More precisely, we want to find Bell inequalities containing only one- and 
two-body correlators. To further simplify the problem we also demand that additionally they obey certain symmetries. From a 
more fundamental point of view, we also want to answer the question of what is the minimal amount of information that it is 
required in order to detect non-locality. 

We shall simplify P in two steps: one is to reduce the order of the correlators (1), and the other is the application of a symmetry 
group G such that the obtained Bell inequalities remain invariant under the permutation by any element of G. By reducing the 
amount of correlators in (1) to those that involve at most K parties, the dimension of P is reduced from (m(d — 1) + 1)" — 1 to 

( 12 ) 

k =1 k ' 

and the number of vertices is kept the same. Each term in the above sum corresponds to correlators of exactly k parties - there 
are (™) ways to choose them. When the k parties are chosen, each of the members can choose one of to measurements, with 
d — 1 independent probabilities (the o'-th one is determined from the fact that the sum of probabilities of all outcomes is equal 
to 1); this amounts to the factor m k (d — l) k . We denote by P k the polytope obtained by not including the correlators of order 
higher than K. 

On the other hand, one can look for Bell inequalities that are invariant under the action of some symmetry group G which 
permutes the parties. The reason to consider this approach that imposes extra constraints is twofold: firstly, even if Vk is simpler 
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than the whole poly tope, its number of vertices is still exponential and its dimension goes like 0(n K ), which can still grow fast 
with n. Adding the permutation symmetry produces a simpler polytope. Secondly, as we will see below, imposing permutation 
symmetry allows for deriving Bell inequalities that can be measured with techniques readily accessible in experiments. 


A. The symmetrized polytope 


Let us consider a group G, which is a subgroup of the group of permutations of n elements, denoted & n . Let us denote by 
Vk the set of all A-body conditional probabilities: 

K 

V K = U • --a>i k \ x ii • ••£**)> 0 < n < ... < i k < n, 0 < dj < d- 1, 0 <Xj Cm}. (13) 

k =1 

G acts on Vk by permuting the indices of the parties i. \ ... fa. Formally, we have defined an action g by: 


g : G x V K —» Vk ^ 

(o’ j Pi 1 ...i k i a i 1 ■ ■ ■ a ik\ x ii ■ ■ ■ x ik)) l— ^ V r 7 (i 1 )...a(i k ) i a ii ■ ■ ■ a i k \ x ii ■ ■ ■ x i k )- 

Every action of a group on a set induces a partition of the set into orbits. Given a probability P £ Vk, we denote its orbit by 

[P] = {g(a, P), a £ G} and the set of orbits by Vk/G. Note that a given probability P can only belong to one orbit. Thus, we 

can express Vk as the disjoint union of all its orbits. 


V K = [J [P]. (15) 

[P]&r K /G 

The sum of all elements in [P] is, by construction, invariant under the action of any element a £ G. Let us then define 
G-invariant probabilities as sums of probabilities that are invariant under the action of G, that is, 

S[p] oc (16) 

Pe[P] 

for all [P] £ Vk/G (there are clearly as many G-invariant probabilities as elements in Vk /G). Accordingly, we define the 
symmetric polytope of /\-body probabilities, denoted P^’, as the image of P under the projection (16). Its dimension is the 
number of G-invariant probabilities obtained from Vk/G. Note that in (16) we have defined the G-invariant probabilities up to 
a proportionality constant. This, however, does not change the shape of the polytope ( i.e. its number of vertices and/or facets and 
their arrangement). Note that P jl is indeed a polytope, because Eq. (16) defines a linear projection. So, P ( / can be completely 
characterized by listing its vertices. However, we are interested in finding the vertices of P^- directly, without first constructing 
the vertices of Vk (of which there is an exponential amount) and then projecting each of them by means of (16). 

As the projection (16) from P k to P^- is a linear map, every vertex of P^ has to be the image of a vertex of P k under this 
projection. The converse is not true in general, because a vertex of Vk can be mapped by (16) to the interior of P^. Since 
every set of deterministic local strategies that give the same value in each ( 16) is a candidate to a vertex, counting the number of 
them will give us an useful upper bound to the number of vertices in P^. To this end, we will make use of Polya’s enumeration 
theorem [62], 

Every deterministic local strategy assigns to each party a list of the predetermined outcomes for all its measurements. So, it 
can be thought of as a function / : X —> Y, where X is the set {(). 1 ,..., n — 1}, which indexes the parties and Y is the 
set of 777,-tuples (y 0 ,... ,y m -i) with 0 < yt < d, in which y,; is the outcome assigned to the z-th measurement. The set of 
all /’s is denoted Y x . Permutation of the parties defines, in a similar fashion as in (15), a partition Y x /G of the set of local 
deterministic strategies: two deterministic local strategies, fi and fa, belong to the same class in Y x /G if and only if there 
exists a permutation a £ G, such that when applied to the parties in fa, we obtain fa. Polya’s enumeration theorem states that 
the number of such classes is given by 


I Y x /G\ = p E l y ! C(<T) ’ < 17 ) 

where c(a) counts the number of disjoint cycles of the permutation a. Because every permutation decomposes into a product of 
disjoint cycles uniquely (up to permuting the cycles), c(a) is well defined. 

Both the dimension and the number of vertices depend on the group G considered, and there is clearly a trade-off between the 
number of elements in G, the number of vertices, and the dimension of P^: the bigger the symmetry group G, the smaller the 
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dimension of P^- and the smaller the number of vertices. 

Example Let us illustrate this result with an example, in which G is the group generated by the shift to the right: 0 i—> 1 i—»• 
... i —} n — 1 i— y 0. In this case, one obtains Bell inequalities which are translationally invariant [63]. One can show that 

\Y x IG\= X ~Y J ^)d mnl \ (18) 

v\n 

where the sum runs over all divisors v of n, denoted v\n, and : ~p(v) is Euler’s totient function, which counts how many integers 
l with 1 < l < v are coprime with v, (i.e., their greatest common divisor is gcd (u, l) = 1). Interestingly, when considering PSf 
in some scenarios, the bound (18) is tight [e.g., (3, 2, 2) and (5, 2, 2)], whereas in some others it is not [e.g. (4, 2,2)]. 


B. The permutationally invariant Bell polytope 

For the scope of this article, however, we shall restrict ourselves to the scenario where the obtained Bell inequalities are 
invariant under any permutation of the parties; i.e., the case where the symmetry group G is & n . Let us first upper bound the 
number of vertices of P®" by explicitly computing (17). In this case, recall that the number of permutations of n elements with 
k disjoint cycles is given by the unsigned Stirling number of the first kind s(n, k), where s(n, k) is defined through the following 
formula [64]: 


E , ,, h / ,\ , (x-\-n —1)1 

s(n, k)x k = x(x + 1) • • • (x + n — 1) = — j- —ryj— 


where a; is a natural number. It then follows that 


\Y X /&u\ 


1 

n! 


E ( rfm ) 

cr£© n 


c(o-) _ J_ 
n! 


J2s(n,k)(d m ) k 

k =1 


1 (n + d m — 1)! 
n\ ( d m - 1)! 


fn + d m — 1\ 
^ d m - 1 ) ’ 


(19) 


( 20 ) 


where the second equality is a consequence of the definition of |s(n, k) |, and the third equality stems from Eq. (19). As we shall 
see in section III D, the bound (20) is not tight in general and it can be refined. 

There is a simpler argument showing that the number of vertices of P®" is upper bounded by (20), which in addition allows 
us to parametrize the vertices of it with d m — 1 integer parameters. The key idea is to observe that a local deterministic strategy 
/ : X —> Y is just a coloring of a hypergraph whose nodes are the elements of X and with the colors being the elements of 
Y. This graph has multiple hyperedges, each one corresponding to a fc-body correlator [cf. Fig. 1], If its nodes are permuted 
according to an element of G, then all the correlators (16) will have the same values. Hence, they will correspond to the same 
vertex of P^-. 



FIG. 1. (color online) The action of a permutation a =* (A C B A) (D i->- D) on a deterministic local strategy. We pick as an 
example the (4, 2, 2) scenario and consider the local deterministic strategy / given by f(A ) = (+, +), f(B) = (+, —), /(C) = (+, +) and 
/(D) = (—, +). Each square represents a party and each circle corresponds to a measurement. We take the convention that if the circle is 
filled in black, its predetermined outcome is +1, whereas if it is filled in red, then its predetermined outcome is —1. All parties having the 
same deterministic local strategy have squares with the same color. Each line corresponds to a two-body correlator; if it is a black line, the 
results are correlated and if it is a red line, they are anti-correlated. Although the value of the single two-body correlators can be changed by 
a (e.g., A o-Bi), observe that the symmetrized correlators do not change their respective value. Both before and after the action of a, one has 
% o] = A 0 + B 0 + Co + D 0 = 2, SjAj] = 2, S[a 0 b 0 ] = 0, Sp 4 0 Bi] = 4, SjAisq = 0. In terms of P, the two graphs above correspond to 
two different vertices, but they are projected to the same point in P®". 
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However, in the case that G = & n , the G-invariant probabilities (16) remain the same for any permutation of the parties. 
Thus, their value depends only on the amount of parties assigned to each color or, equivalently, the amount of parties with the 
same deterministic strategy. There will be as many equivalent strategies as partitions of n in d m (some possibly empty) groups. 
Combinatorically, a simple way to count them is to imagine a set of n elements in a line, then extend it to a set of n + d r " — 1 
elements and pick d m — 1 of them that act as a separator. This will produce a partition of it into d r " groups, some of which may 
be empty. The number of ways to perform this choice is precisely (20). 

From now on, we focus on the case when all the observables are dichotomic ( i.e ., d = 2) and we switch for convenience to 
the expectation-value framework (10), instead of working with probabilities (1). Recall that for d = 2 these two frameworks are 
equivalent. Let us then calculate the dimension of P®", which is the number of elements in P^/©„. Equation (16) directly 
shows that S 1 ^], where Ad is a /c-body expectation value, depends only on the choice of k measurements out of the to, which 
are available. For instance, S\ y \ 2 p, 2 c 0 ] = S[b 2 d 0 e 2 ]- So, one can take as a canonical representative the lexicographically lowest 
element, which in this case would be S[a 0 b 2 C 2 ]’ and simply denote the corresponding correlator as So22- It is then clear that the 
amount of elements in P k/& u is equal to the number of ordered sequences 0 < x\ < X 2 < • • • < x k < m for k = 1,..., K, 
that is. 


K 


\V K /&n\=Y, 


k =1 


k + m — 1 
to — 1 


to + K 
K 


- 1 . 


( 21 ) 


In order to have the simplest symmetric polytope, we will further assume that K = 2, in which case the dimension of P®" is 
just 5. Let us then define the one- and two-body @„-invariant correlators (called simply permutationally invariant correlators) 
we shall be working with: 


n— 1 

S k 0 < fc < 1, (22) 

2=0 


n —1 

S k i := ^ (MfM[ j) ), 0 < k < l < 1, (23) 

i,j =0 

and observe that they are indeed proportional to their counterparts defined in (16). 

Let us now calculate the value of (22) and (23) at every local deterministic strategy. In the case of single-body mean values, 
one can see from Fig. 1 that their value corresponds to the black circles minus the red circles, whereas in the case of the two- 
body correlators one would take the corresponding black edges minus the red ones. It is then natural to work with the following 
variables: 


«/:=ir 1 (+,+)l, bf := |/ _1 (+,-)|, c / :=|/- 1 (-,+)|, df := |/ -1 (-,-)|, (24) 

where / -1 is the preimage of /. Notice that a / counts how many of the parties have predetermined outcomes Ado = +. Af -| = + 
for the deterministic local strategy /, etc. In the example of Fig. 1, one has aj = 2, bf = 1, c/ = 1 and df = 0. Observe 
that, for all /, one always has that «/ + bf + Cf + df = n because there are no more possibilities. Interestingly, this is all the 
information needed to compute the value of (22) and (23) at any deterministic local strategy. Indeed, one directly sees that, for 
a given /, <5>o = a/ + bf — Cf — df and Si = af — bf + Cf — df. Moreover, it follows from Eq. (11) that the value of the 
two-body permutationally invariant correlators (23) for any / is 


n— 1 

S kl = S k St - Y, (25) 

2=0 

where for k = l the subtracted sum amounts to n , while for k ^ Z it is given by 

n— 1 

2 := Y( M o )M i } ) =af~ b f^Cf+ df- (26) 

2=0 


Finally, we have that <Soo = (<^o) 2 — <5oi = SqS± — Z and Su = (<Si) 2 


— n. 
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C. Permutationally invariant Bell inequalities with two-body correlators 


The general form of Bell inequalities constraining the polytope P 2 in the scenario (n, 2, 2) is 


n— 1 

2=0 


+ E E S ij (M^M[ j) }+ E ^(M^M^)>-p c , 

0 <i<j<n 0 <i^j<n 0 <i<j<n 

(27) 


for some «i, /3*, 7 ij, Sij,£ij £ R and /3 C £ R, which is the so-called classical bound defined as Pc = — min^ e p 2 J 2 . There are 
2n 2 degrees of freedom, which can be read from expression (12) for m = K = d = 2. 

Then, the general form of Bell inequalities bounding the permutationally invariant polytope P®" simplifies to 


Pc + OtS 0 + PSi + T^oo + &S01 + — (Sn — 


(28) 


where a , /3, 7 , <5, £ £ R and p c £ R is the corresponding classical bound. In other words, by symmetrizing over 6 n , we have 
found a way to obtain Bell inequalities of the form (27) with all the coefficients being equal: ct* = a, Pi = /?, 7 ^ = 7 , S t j = 6, 

and £j j = e. 


D. Characterizing the vertices of the symmetric polytope 


The symmetrized 2-body polytope P®" in the (n. 2, 2) scenario is embedded in a 5-dimensional space with coordinates 
(So, (Si, (Soo, (S 01 and (Sn. However, if a given set of correlations corresponds to a local deterministic strategy /, only 4 parameters 
are needed in order to completely describe it: either (a, b, c, d) or (n, S±,So , Z). To simplify the notation, we have omitted the 
subindex /. Interestingly, these two sets of variables are related via the following (unnormalized) orthogonal transformation 




where H stands for the 2x2 Hadamard matrix 





(29) 


This relation can be generalized to the (n, m, d ) scenario with if-body correlators, both in the framework of expectation values 
( 10 ) and in the framework of probabilities ( 1 ), but this will be studied elsewhere. 

In order to characterize all the vertices of P® n , let us introduce the following set of all the four-tuples in the set of non-negative 
integers that sum up exactly to n: 

T n := {(a, b, c, d) £ Z 4 : a, b, c, d > 0, a + b + c + d = n,}. (30) 


Geometrically, T„ is the set of points of a simplex with integer coordinates. In the (n, 2, 2) case, it can be viewed as the set of 
points of a tetrahedron with integer coordinates. Then, the polytope P f" can be defined as the convex hull of y?(T n ), where <p 
is the following map: 

(p: T n > R 2 " 

(a,b,c,d) 1 —> ((Sq, (Si, (Sqo, (Sqi, (Sn) 


Importantly, one can discard most of the points in T„, as they get mapped onto the interior of P®". In fact, as we prove in the 
following theorem, the vertices of P®” are completely characterized by those four-tuples that belong to the boundary of T n . 
Let us denote by Ext(P®") the set of vertices of P®" and let us define 

9T„ := {(a, 6 , c, d) £ TLp : a,b,c,d> 0, a + b + c + d = n, abed = 0}, (32) 


which is the set of points of T n with at least one coordinate equal to 0 . 

Theorem 1. For all p = (a, 6 , c, d) £ T„, the following equivalence holds: 

ip{p) £ Ext(Pf"). 


p £ <9T„ 


(33) 
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A proof of this can be found in [52]. For completeness, we have included it also in Appendix A. Theorem 1 enables us to 
construct P®" by taking the convex hull of ip(cPT n ) instead of the convex hull of ip(T n ). Hence, the total number of vertices is 

|Ext(P®")| = |0T n | = |T n | - |T n \ 8T n \ = |T n | - |T„_ 4 | = ^ ^ = 2(n 2 + 1). (34) 

Hence, we have improved the bound from 0(n 3 ) in (20) to 0(n' 2 ) in (34). The bound (34) is now exact due to Theorem 1 and 
the fact that ip is invertible, so any vertex of P®" can be generated from a unique tuple (a, 6, c, d) £ <9T„. We may interpret this 
in terms of the graph introduced in Figure 1 ; it means that a coloring in which all the different colors appear cannot correspond 
to a vertex of P®”. Thus, the coordinates of all the vertices of P f" are of the form 


Sq = a + b — c — d 

(35) 

Si = a ~ b + c — d 

(36) 

5oo = (So) 2 — n 

(37) 

S 01 = 6>oiSi — Z 

(38) 

Sn = (Si)- — n, 

(39) 


for (a, b, c, d) £ dT n . 

IV. CLASSES OF BELL INEQUALITIES 

Here we show in detail how the Bell inequalities presented in [52] were found. Later we will also characterize its tightness. 

In order to search for particular classes of Bell inequalities, we will use the parametrization we have derived in (35-39). It 
follows from Theorem 1 that P®" = conv(<^(<9T n )), where conv stands for convex hull. However, if one relaxes the condition 
on the numbers (a, 6 , c, d ) from being integer to being real, one obtains a new convex object — not a polytope, but easier to 
characterize. Let us consider the following sets: 

T„ = {(a, b,c,d) £ R 4 : a, b, c, d > 0, a + b + c + d = n} (40) 

and 

dT n = {(a, b, c, d) £ R 4 : a,b,c,d> 0, a + b + c + d = n, abed = 0}. (41) 

We shall define P® 71 = conv(<p(T„)). Since the proof of Theorem 1 also applies to the continuous case, conv(y>(T„)) = 
conv(<y 9 ( 9 T„)). Because <9T„ C dT n (or T„ C T n ) one has then P® 7 * C P,® 7 \ Consequently, a Bell inequality valid for P®" 
is also a Bell inequality valid for P®'". We expect that finding Bell inequalities of the set P®" is an easier task than finding the 
ones corresponding to P®", at the price of having inequalities not as optimal as the ones of <9T„ C dT n . 

In passing let us notice that the sets <p(d T„) and are semialgebraic, since they are defined through polynomial equali¬ 

ties and inequalities. The characterization of convex hulls of semialgebraic sets is a well studied subject [65]. Although the exact 
characterization is an NP-hard problem [65], there exist efficient approximations with semi-definite programming techniques in 
terms of the so-called theta bodies. These techniques can be applied to any (n, in. d) scenario with A'-body correlators when 
the symmetry group G is S„ and they will be studied elsewhere. These approximations are in the spirit of the Navascues- 
Pironio-Acin (NPA) hierarchy [57], which bounds the set of quantum correlations Q appearing in ( 8 ) from above by providing 
a hierarchy of sets QC ... C Q 2 C Qi+ab C Q i . Each level of the hierarchy better approximates Q and Q/ ;: asymptotically 
converges to Q as k —> oo. When the NPA hierarchy certifies that a set of correlations is outside Q/.. for some k, then these 
correlations cannot be realized with quantum resources. In our case, a characterization of P® " through the so-called theta bodies 
produces also a hierarchy of sets P®" C ... C 0 2 C 0 4 that in this case certifies that a set of correlations which is outside of 
(-)for some k cannot be simulated through shared randomness and it must be necessarily nonlocal. 

The approach that we shall follow here in order to derive a class of Bell inequalities is based on the following: Eqs. (35-39) are 
polynomials of at most degree 2 in the variables a, 6 , c and d. They define a 3-dimensional manifold (because of the constraint 
a + b + c + d = n) in a 5-dimensional space. By finding a tangent hyperplane to it, i.e., a plane of the following form: 

a5o + PS\ + t^oo + SSo! + 2^ 1:l + Pc = 0, (42) 

and then proving that the LHS of (42) is positive on T„, we shall obtain the coefficients a , ;3, 7 , <5, e and the classical bound (3 C 
that define a valid Bell inequality. Theorem 1 also hints us to look for facets of T n , that is, points (a, 6 , c, d) with coordinates 
satisfying abed = 0. We aim at obtaining a Bell inequality which is tangent on vertices of P®" represented by four-tuples 
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belonging to a single facet of T n , for example the one defined by a = 0. Let us, however, remark that not all Bell inequalities 
need to come from the same facet, as a = 0 implies abed = 0 but the converse does not hold, so the class of inequalities we 
shall derive will display this feature - this is discussed in detail in the next sections. 

To be more rigorous, we consider the following Lagrangian function 

£ = aSo + /3S\ + < ^ 01 3” AbAcAd — n)-\- p,a , (43) 

and we wish to find its minimum — /3 C = min aj f, jC ,d,A,MeiFt £■ The necessary condition for an extremal value of £ to exist in the 
point x* = (a* , b*, c *, d* , A* , p,*) is that the following system of differential equations 

d£ 

— =0 (x = a,b,c,d, X,p) (44) 

is satisfied in x*. Since all the permutationally invariant mean values appearing in Eq. (43) are at most quadratic in the variables 
a, b , c, d, one finds that (44) can be rewritten as a simple system of linear equations 


/ £ + c -c -tf 1 1\ 


( a* \ 


/ —a — /3 + C +/2 \ 

c £_ -c 1 0 


b* 


-a + (3 + £_/2 

-c Z- CIO 


c* 


a—/? + £_/2 

-c+ -c C C+ 1 0 


d* 


cr + /3 + C +/2 

1 1 1 100 


A* 


0 

V 1 0 0 000/ 


V M* / 


V 0 ) 


(45) 


where := 7 ± 2 5 + e and £ := 7 — e. 

For a generic choice of and £, the matrix from (45) is non-singular, leading to a unique solution (a*, b*, c *. d*), i.e., a 
Bell inequality which is tangent to P®" only at a single point. This is what one would expect from a generic choice of the 
coefficients a, /?, 7 , <5 and e. And, although such solutions would lead us to proper Bell inequalities, we would like to find Bell 
inequalities that are tangent to P®" at more points, and therefore stronger in detecting nonlocality. For this purpose, we can 
impose some constraints on the coefficients 7 , S and e appearing in Eq. (45). Precisely, the number of solutions of Eq. (45) is 
given by Rouche-Frobenius theorem, which says that a linear system of equations Ab = c is compatible if, and only if, the rank 
of A is the same as the rank of the extended matrix A\c; the solution is unique if, and only if, det A f 0 . As we do not want a 
unique solution, the first condition is that det A = 0, which when applied to Eq. (45) is 


d 2 — 7£ = 0. 


(46) 


With this condition the rank of A is not greater than 5, and it is exactly 5 if we further impose 7 6 s f 0, which we can assume to 
be always fulfilled. Otherwise, Eq. (46) would just produce trivial Bell inequalities. In order to guarantee a solution of (45), we 
have to ensure that also rank(A|c) = 5. This is done by imposing that the 5x5 minors of A\c vanish, which is equivalent to the 
following condition 


6(/3 A S) = s(a A 6). 


(47) 


In this case, the classical bound reads 


/3 C = —£(a* 


b*,c*,d*,\*,p*) 


(/3 + S) 2 A n(5 — e ) 2 
2s 


(48) 


If the minimum is global, we have obtained a Bell inequality which is tangent to P.®” and valid for P®", although in general 
not tangent to the latter. By repeating the same argument, but over the facets b = 0, c = 0 and d = 0 we would get similar 
expressions as in (47) and (48), with swaps of parameters and changes of signs. In the following theorem [52], we generalize all 
of them and we do a small correction to (48) ensuring that the Bell inequality is tangent to P® n , and not to P®". We also prove 
that the minimum is global, i.e., the Bell inequality is valid. 

Theorem 2. For any cr £ {— 1,1} and x,y,^L £ IN such that y has opposite parity to x (y) if n is even (odd), define the 
parameters of (28) as 

a± = x[ap± (x + y)\, (3 = p,y, j = x 2 , 6 = crxy, e = y 2 . (49) 


Then the classical bound of the resulting Bell inequality, for which it is tangent to P 


&n 
2 ’ 


is 


/3c 


1 

2 


[n(x + y) 2 + (ap, ± x) 2 


!]■ 


(50) 
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Theorem 2 shows that the few-parameter class of Bell inequalities defined by Eqs. (49, 50), which can be explicitly written as 

2 ^ 1 

: = x \ a B ± ( x + y)) s < 0 + M yS-L T yS 0 o + oxySoi T yiSn > —-[n(x + y ) 2 + ( 07 / ± x) - 1], (51) 

is tangent to P ®in at least one vertex, but it does not guarantee that it is tight (that is, defines a facet of P®") and as such 
optimal. Notice that in this case proving that a Bell inequality defines a facet of P®’* requires showing that the set of vertices of 
P® n saturating it spans a four-dimensional affine space. In the next theorem we are going to characterize the vertices of P®" 
which saturate the Bell inequality (51). To this end, it will be useful to rename the faces of the tetrahedron T„ according to the 
choice of a and ± in order to obtain a more general result: 


a 

+ 

+ 

- 

- 

± 

+ 

- 

+ 

- 

r 

b 

c 

a 

d 

s 

a 

d 

b 

c 

t 

d 

a 

c 

b 

u 

c 

b 

d 

a 


Now, Eq. (29) is generalized to 


( n ^ 


(i 



( r \ 

( 

Si 

= 2 

To¬ 



S 

= 

So 


il 



t 





-a ) 


\u) 

V 


i 

1 

l 

l 

To¬ 

b 

-H 

-H 

b 

-H 

il 

±1 

Tl 

Tl 

—a 

G 

G 

— G 


\ 

( r 


s 


t 

) 

\u) 


(52) 


(53) 


where a and ± are independent signs and // is the 2x2 Hadamard matrix defined in Eq. (30). Having this, we can now state 
the following theorem. 


Theorem 3. Let us assume that x and y in the Bell inequality (51) are coprimes and let us define the following quantities: 


K r (r ) := [±n(y — x) + crp ± x + t\/ 2, 
to(r) := ±y~ 1 K r (r) mod x, 
uo(t) := [nx ± K r (r) - (x + y)t 0 (r)]/x, 
So{r) := [TA' r (r) + yt 0 (r)]/x, 


where r £ {—1,1}- If a vertex of P®" saturates (51 ), then it is of the form <p([r, s, t , u]), where 

[r, s, t, u ] = [0, s 0 (t), t 0 (r), u 0 (r)] + fc[0, x, y, -(x + y)], k £ TL. 
Furthermore, the number of vertices of P®" saturating (51 ) is given by 

Uo(r) 


Ns := ^ max < 0, 


r=±l 


x + y. 


— max < 0 


-«o(t) 


+ 1 


(54) 


(55) 


(56) 


Clearly, a necessary condition for the inequality (51) to be tight is that Ns > 5 as 5 is the minimal number of points to generate 
a 4-dimensional affine subspace. We have checked numerically that this condition is also sufficient for all the instances of n for 
which we could compute all facets of P®", which is n < 33. To prove sufficiency, let us make the following observations: The 
proof of Theorem 3 indicates that the vertices saturating (51) form a zig-zag pattern that alternates between r = 1 and r = — 1. 
For a given t, more than three different points are always linearly dependent. Hence, one must pick three points for r and two for 
—r. Short algebra shows then that the matrix formed by the coordinates of such points (in particular, if chosen with consecutive 
k) has maximal rank. 


V. QUANTUM VIOUATION 

In this section we analyze the quantum violation of permutationally invariant 2-body Bell inequalities for the (n, 2, 2) scenario; 
i.e., Bell inequalities of the form (28). If the left-hand side of (28) is smaller than 0 for some //, then the inequality is violated and 
the nonlocality of this p is detected. To show that there is a quantum violation of such an inequality, one must find a quantum 
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state (taken from a Hilbert space of some dimension D) and POVM elements such that the probabilities obtained through (6) 
violate (28). 

In order to find out what is the best performance of a Bell inequality, we are interested in finding the states and measurements 
giving the maximal quantum violation of such a Bell inequality, however, the dimension D of the corresponding state depends 
in general on the scenario under consideration. Luckily, in the (n, 2, 2) scenario, D = 2 is enough, because in order to find the 
maximal quantum violation of a Bell inequality with two dichotomic measurements per party, it suffices to consider qubits and 
traceless real observables [66, 67]. 

Let us define by a = [a x ,a y , a z \ the usual vector of the Pauli matrices and let us denote by c?W the same vector acting on the 
?'th site of an n-partite system. Any traceless qubit observable with eigenvalues ±1 can be written as n ■ a with n £ R 3 such 
that |n| = 1. Then, any real traceless single-qubit observable Mix] with eigenvalues ±1 can be parametrized as 

Mg = co s(9g)ag + sin (9g)ag, (57) 

where 9g £ [0, 27 r). To simplify the notation, we shall denote Og as p>i when Xi = 0 or as 9i when x z = 1. 

The so-called Bell operator is the quantum observable corresponding to the left hand side of (28), where the observables (57) 
are used in the definitions of (22, 23). We shall denote it B({pi, 9i}). The maximal quantum violation is then given by the 
minimal eigenvalue value of B({pi, 0,}), when varying the angles ipi and (/,, and the corresponding eigenvector is the quantum 
state maximally violating the corresponding Bell inequality (28). 

If a Bell inequality displays some symmetry with respect to an exchange of parties (i.e., it consists only of G-invariant 
correlators as defined in (16) for some group G), the pure state | ip) for which the maximal quantum violation is achieved does 
not need to display the same symmetry (see [63]), nor do the observables be the same at each site (i.e., Mlg be independent of 
i). However, one can always consider a mixed state p and observables (at the expense of increasing the local dimension D) such 
that p is invariant with respect to any permutation in G and the extended observables are the same at each site. However, in the 
case we are considering, G = ©„, numerical results for a small number of parties suggest that the maximal quantum violation 
is already achieved with a permutationally invariant pure qubit state and the same pair of measurements at each site. This is 
very convenient for our analysis because in the latter case the Bell operator B(p, 9) depends on only two angles p = p, and 
9 = 9i, and it is permutationally invariant and therefore can be block-diagonalized. Finally, as the Bell inequality (28) contains 
only 2-body correlators, we are able to prove that in a certain basis, B(p, 9) is pentadiagonal. All this implies that its minimal 
eigenvalue over p, 9 can be found numerically for a large number of parties in an efficient manner due to sparsity of the matrix 
and the fact that we have reduced the optimization to a single variable problem. 


A. Block-diagonalization of the permutationally invariant Bell opeator 


When the same pair of measurements is taken at each site, the Bell operator is invariant under any permutation of the parties. 
It is a well-known result that under these circumstances it can be block diagonalized. The reason for that lies in the result known 
as Schur-Weyl duality [68]. Consider the permutation operator V a acting on an n-qubit Hilbert space: 

1G . |io * * * i n _i) i y Iv-ifo) * * * i(j— !(n— i)}^ (58) 

where a £ & n . V a is a unitary representation of a £ & n acting on the n qubit Hilbert space. Consider also the tensor product 
representation W u = £/® n of the special unitary group. The Schur-Weyl duality states that the Hilbert space splits into blocks, 
on which these two representations commute (in representation theory terms, the action of these two groups is dual). This block 
splitting is the following: 

n/2 

((D 2 )®" = (J) Hj ® /Cj, (59) 

J=J 0 

where 2 J 0 = n mod 2, the Hilbert spaces 77, / are of dimension 2 J + 1 and the dimension of ICj is 1 if J = n/2 and 


dim/Cj = ]n/2 — j) ~ (ra/2 — J — l) <60) 

otherwise. 

A permutationally invariant state p is such that for all a £ ©„, p = V a pVJ. In the basis given by (59), p has the following 
form: 


P = 


n/2 

cr\ pj 

jSo dim/C ' 


PJ & ldim/Cj ? 


(61) 
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where pj is a probability distribution and pj are density matrices. Hence, p, which is a 2" x 2 ra density matrix, is completely 
characterized by the [n/ 2J blocks pj, the biggest one having size (n + 1) x (n+ 1), corresponding to the symmetric subspace 
id' 2 )-". 

The blocks pj can be obtained by projecting the Bell operator onto the following basis [69]: 

{\D k 2J )®\rr m/2 }iu, ( 62 ) 

where m = (n — 2 J) is an even number, | D’f) are the so-called symmetric Dicke states (of n qubits and k excitations) which 
are symmetric superpositions of k qubits in the state |1) and n — k qubits in the state |0): 

/ \ — 1/2 

1 D n) : =r k ) E ^(io)® n_fe ii)® fe )> (63) 

' ' ae&n 


and I'i/j ) is the singlet state given by 


IV’ ) 


| 01 )-| 10 ) 

Vi 


(64) 


In the following theorem, we give the analytic form of each block, when the same set of real traceless observables is taken at 
each site: 


Mq' 1 = cos pa^ + smpa^\ = cos 9a^ + smOa£\ 


(65) 


Theorem 4. The J-th Block Bj(p,9) of the Bell operator corresponding to a 2-body symmetric Bell inequality (28), with 
measurements given by (65), has elements ( Bj(p, 9 ))[' (k-th row, l-th column) given by 


(Bj(<p, 0 ))f = dkSk,i + UkSk,i -i + uiSk-i,i + VkSk,i-2 + viSk-2,1, 


( 66 ) 


where 0 < k,l,< 2 J, dk (0 < k < 2J) are the diagonal elements, Uk (0 < k < 2J — 1) correspond to the elements of the 
upper (lower) diagonal and Vk (0 < k < 2 J — 2) stand for the elements of the second upper (lower) diagonal and 5k,i is the 
Kronecker delta function (5k,i = 1 4=> k = l, otherwise 5k,i = 0). The coefficients dk , Uk , Vk are explicitly given by 


d k ■= P c + (2 J - 2 k)A + [(2 J - 2k) 2 - n\B/ 2 + [2fc(2J - k) - m]C/ 2, (67) 

Uk ■= [A’ + (2J - 1 - 2k)D] ^(2J ~k)(k + 1), (68) 

v k := C\J (2J - k)(2J -k- 1)(A: + l)(fc + 2)/2, (69) 

where the parameters A, A' : B , C , D depend only on the Bell Inequality coefficients and the measurements’ angles: 

A \= a cos tp + (3 cos 0, (70) 

A' := asin^J + /3sin0, (71) 

B := 7 cos 2 p + 25 cos p cos 0 + e cos 2 0, (72) 

C := 7 sin 2 <p + 25 sin p sin 9 + e sin 2 6 , (73) 

D := 7 cos p sin p + 5 cos p sin 0 + 5 cos 9 sin p + e cos 9 sin 0. (74) 


Observe that (66) ensures that every block of the Bell operator is pentadiagonal, and this comes from the fact that we are 
considering a symmetric Bell inequality with at most 2-body correlators. From the proof of Theorem 12 in Appendix B, it 
is clear that if the Bell inequality involved 3-body correlators, it would be heptadiagonal; if it involved 4-body correlators, 
enneadiagonal, and so on. This argument also shows that nonlocality of some of the highly entangled multipartite states cannot 
be revealed through Bell inequalities involving only few-body correlations: if we consider the GHZ state (l/v / 2)(|0)®" + |l)®”) 
which is supported on the last block of (61) corresponding to J = n/2, the density matrix of the state has only four terms, two 
of them lie in the diagonal and the other two are the coherences corresponding to |H°)(D™|, |H”)(i9°|, which can only be 
reached through a full-body correlator. Hence, for a symmetric 2-input 2-output inequality, which is not full-body, the GHZ 
state is indistinguishable from the classical mixture (l/2)(|0)(0|® n + |l)(l|® n ). However, there are other states, which show 
nonlocality in a robust way, some of them even experimentally realizable, as we shall discuss it in further sections. 

The form (66) also has important numerical implications, as the Bell operator can be easily stored in a sparse matrix. This 
allowed us to show nonlocality of genuinely multipartite entangled states of more than 10 4 qubits [52]. Recall that, in order to 
find the largest quantum violation with these settings (i.e., quantum state and measurements) of a particular Bell inequality, we 
have to vary the angles p, 9 in (66) and look for the smallest eigenvalue of the [n/2j possible blocks. Actually, as it is proven in 
the following fact, the spectrum of the Bell operator depends only on the difference between p and 9. 
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Theorem 5. Let B(p, 9) be the permutationally invariant Bell operator corresponding to the measurements given in Eq. (65). 
Then, for any cfE, B(p + c, 9 + c) = \U(c)\® n B(<p, 0)[f7(c)']®” with U(c) being a unitary operation given by 


U(c) := 


cos(c/2) — sin(c/2) 
sin(c/2) cos(c/2) 


(75) 


Theorem 5 implies that the spectrum of B(p,9) depends only on <p — 9. This fact will later be used to simplify our analysis 
of the maximal quantum violation of the corresponding Bell inequalities. Another simplification comes from (66) and the form 
of the constant C given in Eq. (73). Concretely, there is a convenient choice of p and 9 for which C = 0, which makes the Bell 
operator B(p, 6) tridiagonal. Such form is given in the following corollary: 

Theorem 6. Let k = 9 — ip. Then the Bell operator B(p, 6) is tridiagonal when 


9+ = arctan 


7smre 
7 cos k + 6 ± \J5 2 — ye 


B. Analytical class of states 


Here we present an analytical class of states for which the Bell inequalities from the family (51) are violated. We also show 
that this analytical class of states converges to the maximal quantum violation when n becomes large, providing strong numerical 
evidence for that. Such a state is a Gaussian superposition of Dicke states with variance growing as 0(y/n) and centered at the 
minimum of the function dk defined in Eq. (67). We have found that for this class the maximal quantum violation is achieved 
with a state from the symmetric block (J = n/2), although other blocks also lead to quantum violation. In other scenarios, such 
as the ( n , 3,2) scenario, we have also found some Bell inequalities maximally violated with states supported on the lowest J 
blocks. 

Theorem 7. Let \ ip n ) = \^n) b e a symmetric n-qubit state with the following coefficients in the Dicke basis: 


/ ( n ) e 
n ■= 


-(fc-M) 2 /4 a 


where p := n/2 + A/(2B — C ) and e 1 /2n -C tx <C n. Moreover, let 9(p) be the solution ofEq. (6) guaranteeing that (7 = 0. 
Then, 


(^n\B n /2(V,0{v))\^n) 


'Pc 

-?+e- 1/8CT 

(a' - 45)1 

n + 

n 

2 

( b ) J 



2Ba -— + e~ 1/8a ( A' 
2 B 



+ o(an 1 ). 


Note that the choice of p is motivated by the following fact that dk is a quadratic function in k with minimum at k = p. 
Hence, if B n / 2 {p, 9) were diagonal, where [•] is the rounding function, would be its minimal eigenvalue. However, since 
this value is always positive for the family of inequalities described in Theorem 2, we must make the off-diagonal elements play 
some role in order to achieve quantum violation. 


Example. Let us illustrate Theorem 7 with an example. Consider the Bell inequality presented in Ref. [52], which is a particular 
case of the class (51) with the choice of parameters x = y = 1, u = ± = — 1 and p = 0, which reads 

— 2Sq + 2^00 ~ <5>oi + 2'^ 11 + > 0. (76) 

We have that C = (siny; — sind) 2 . Hence, the solution of C = 0 which makes Ado ±Adi is to take p = n — 9. Then, the 
values of A , A', B , C, D for a = —2, /3 = 0,7 = 1, 5 = —1, e = 1 and ip = n — 9 are 


A = 2 cos (9, A’ = — 2sin0, H = 4cos 2 0, 

(7 = 0 D = 0. 

Therefore, the form of B n / 2 ^: — 9 , 9) is then given by 

dk = 2n(l + cos 9 + (n — 1) cos 2 9) — [4 cos 0(1 + 2n cos 9)]k + 8 k 2 cos 2 9 , 
Uk = —2 sin Q-\J (n — k)(k + 1), 
v k = 0. 
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Hence, the class of states proposed has parameters 


n 1 

Mn = 7T -< 7)) 

2 4 cos 0 

and <T n £ Q(y/n) such that it fulfils the requirements for the validity of the approximations of Theorem 7. The asymptotic 
expectation value is then given by Theorem 7: 

(<Pn\B n/2 (n - 6,8)\(p n ) ~ ~ ^ +e~ 1/8 ‘ 7 A , ' S j n+ (^Ba - ^ + e~ 1/8<T A'^j To^n" 1 ). 

To find the maximal value of the above expression let us approximate its leading term by 

— ~^ + e- 1/8 M' + A' = 2(1 - cos 2 d - sin 0), 

n 2 n 2 

and notice that the latter attains its maximum for 6* £ {7r/6, 57 t/6}. 

As a result, the quantum violation of inequality (76) by our analytical class of states | ip n ) divided by its classical bound reads 
(let us take a = yfn for simplicity) 

^(v?n|£n /2 (V6,57r/6)|<Pn) - + 3n _1/2 - ^n _1 + o(n~ 3/2 ). 

Hence, the quantum violation exists for any n and it tends to —1/4 when compared to the classical bound as n —► oo. 


C. Accuracy of analytical results 

Let us now compare the analytical results provided above to the maximal quantum violation of the Bell inequality (76) 
obtained numerically. For this purpose, let us denote as It/™ 111 ) the state maximally violating (76). As we have discussed, it is a 
superposition of Dicke states, and thus can be written as 


K“) = 


E4"> 

0 


D k ) 

^nh 


(77) 


(n) (n) 

where c k are real coefficients, because the Bell operator is taken to be a real symmetric matrix. Figure 2 shows that c, display 
an interesting regularity when plotted versus k\ this corresponds to the case H„/ 2 (0 ,0 — <p), where 0 - ‘-p has been optimized 

numerically. However, if ip is such that C = 0, then the form of the cl// now denoted a k 'K is displayed in Figure 3 and 
compared with the form given in Theorem 7, where we can see a good agreement. The relative difference between the numerical 
and the analytical state giving maximal quantum violation is compared in Figure 4. Observe that Theorem 5 ensures that both 

(n) (n) 

c k and a k lead to the same quantum violation; one just has to rotate the measurements accordingly. 


D. The 2-body reduced state 

Using the fact that, when we trace out a number of subsystems of a quantum state supported on the symmetric space, the 
reduced density matrix of the state is the same regardless of which systems are traced out, we are able to give an equivalent way 
of computing the maximal quantum violation of a permutationally invariant Bell inequality (28) provided all parties perform 
the same pair of measurements. Indeed, if we denote by p 2 the two-body reduced state of | ip), where | ip) maximally violates 
inequality (28) with measurements A4o and A4i, then we only need to multiply the expectation values (Aik ® Aii) P2 or 
(A4fc <8> l 2 ) Pa by the times they appear in Su or in Sk■ By doing so, the quantum expectation value of (28) simplifies to 

(ip\B\ip) = fic + n (a(Mo <8 l2) P2 + P(Mi 8 1 2 ) P2 ) + n(n - 1) Q(A4 0 8 M Q ) P2 + 6(M 0 8 A4i) P2 + 8 A4i) P2 ) 

(78) 

And if the value of (78) is negative, nonlocality is revealed. 

Let us first show how to compute this reduced state: 
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N = 10 3 



N = 10 4 



k 


FIG. 2. (Color online) Distributions of c^ versus k for different values of n = 10,10 2 ,10 3 ,10 4 . are the coefficients of the eigenstate 
corresponding to the minimal eigenvalue of B n / 2 ( 0 , 6 — ip), for the optimal 6 — ip. The black dots correspond to the values of c'f * and the 
blue line connecting them is added to better show how they behave with k. 


Theorem 8. Let p be the density matrix of an n-qubit state supported on the symmetric space. Let p$ be the representation of 
p in the symmetric space. Let us denote by Tr„_ c j(p) the density matrix of p after tracing out n — d subsystems. This reduced 
density matrix is given by 


(Tr n - d (p)) 


Li-JUL |)UV,)’ 


(79) 


where 0 < i 7 , j ; < 2 d , i' = iq .. .id- 1 and]' = jo .. -jd-i ore represented in binary and |i'|, |j'| is their number of ones in this 
binary representation. Columns are indexed by and rows are indexed by ]'. 

Let us observe that the factor appearing in Eq. (79), i.e., 


f(n, k, d, |i'|, |j'|) = 


n — d 
k 


n 

k+ |j'| 


” 1/2 


(80) 


is symmetric in |i'| and |j'| and it fulfills the geometric mean property 


f(n, k, d, |i'|, |j'|) = y/f(n,k,d, |i'|, \i'\)f(n,k,d, |j'|, |j'|). 


(81) 


When exploring large values of n. from the point of view of numerics, it is important to have in mind that the combinatorial 
numbers appearing in / grow very rapidly, affecting the formula given in Eq. (79) if they are not handled carefully. In particular, 
for n > 650 they can take values greater than 10 308 , which is the storage limit for 64-bits floating point arithmetic. It is therefore 
important to work with a simplified version of / and here we present the cases for d = 1,2, which are the ones we are going to 
use in the present context. 


n — k 


k + 1 


f(n,k, 1,0,0) 


n 


f{n,k, 1,1,1) 


n 


(82) 
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FIG. 3. (Color online) The black points represent values of a J/ 1 ' for n = 10,10 2 ,10 s , 10 4 . The cyan line is added to mark better their 
behavior with k. The red lines present square roots of the Gaussian distributions of Theorem 7 with means p n and variances <r„ chosen so 
that the distributions fit the points. Their explicit values are p n = n/2 + (1/4 cos 6„), where 6ho = 2.6672, # 10 2 = 2.6334, # 10 3 = 2.6231, 
9 10 4 = 2.6 1 80, and ctio = 0.4049, <j 10 2 = 1.3935, cr 10 3 = 4.5 109, and <r 10 4 = 14.379, respectively. Noticeably, already for n = 10 the red 
line matches the points very well. 


f(n, k, 2,0,0) = 


. — k n — k — 1 
n n — 1 


f(n,k, 2,1,1) = 


n—k—lk + 1 
n n — 1 


f(n,k, 2,1,1) = 


k -f- 1 k -T 2 
n n — 1 


(83) 


Using Eq. (81) we can find the rest of the values. 

We are now ready to show what is the form of p 2 for the state given in Theorem 7. 


Theorem 9. Let \tp) be the state introduced in Theorem 7. Its two-body reduced state p 2 can be well approximated for large n 
as 


/ 

— (1 + 2c) 

—c 

1 1 

\ 

« 2 i+x+r+f 

—c 

—c 

0 

0 

0 c 

0 c 

+ o(n) 

V 

\ 1 

c 

c (2c — 1) ) 

/ 


(84) 


where c := p — n/2 and |+) := (|0) + \l))/\/2. 

It is worth noting that, by monogamy of entanglement, p 2 tends to a separable state, as one can see directly from (84). When 
this state is used in the class of Bell inequalities presented in Theorem 2, we observe that the quantum violation relative to the 
classical bound increases [52], i.e., its quantum violation increases linearly with n, whereas Theorem 9 ensures that the two-body 
reduced states are almost separable. The reason for that is that the coefficients of the inequality in Theorem 2 are picked such 
that the state (|0) + \\))/s/2 has expectation value exactly zero on them, so the quantum violation is completely defined by the 
second term in Eq. (84). 
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FIG. 4. (Color online) Maximal quantum violations of inequality (76) and obtained from diagonalization of the symmetric block B n / 2 (blue 
crosses) are compared to violations by the states defined in Theorem 7 (red crosses). All means and variances in the Gaussian distributions 
were chosen so that |7/;™ ln ) maximally violate (76). The inset contains the relative difference between numerical and analytical violations. 
The points in the inset are the actual values of n for which we performed the calculations and the cyan line is added just to help visualize the 
tendency. Clearly, already for n = 10, the violations are almost the same. 


E. Experimental realization of permutationally invariant Bell inequalities with 2-body correlators 

Let us finally notice that when the two-body Bell inequality is permutationally invariant and the observers perform the same 
pair of measurements at each site (65), the one-body and the two-body permutationally invariant expectation values Sk (22) and 
Su (23) can be determined by measuring collective spin operators. Such operators, corresponding to spin-1/2 particles, are 
given by 


S r 


n— 1 

; = E' 


i=0 


(85) 


where a £ {x,y,z}. It is convenient to consider the vector S := [,S' X , S y , S z \. The value of the permutationally invariant 
correlators (hence, the value of any Bell inequality built from them), can be readily measured by taking measurements in any 
direction n in the Bloch sphere, where n is a unit vector, by means of n ■ S. Observe that, if is aligned along the direction 
rhfc for all i, meaning that = rhfc • o, then Sk is given by 


S k = 2(rhfc • S). 


( 86 ) 
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Similarly, the value of S k i can be expressed in terms of the different directions mo, mi of the measurements as 


Ski 


in. + ih;) 




n(m k ■ i hi). 


(87) 


In particular, when k = l, we have 



( 88 ) 


The way to estimate such collective measurements was already proposed in [78] and it was implemented [106] using spin 
polarization spectroscopy [110]. In section VII we discuss its implementation with other technologies. Moreover, motivated 
by the recent experiment aiming at detection of entanglement in the Dicke states [56], in appendix C we study how certain 
experimental errors affect the quantum violation of our Bell inequalities. 


VI. INEQUALITIES DETECTING DICKE STATES 


Symmetric Dicke states | D%) [70] are of great importance in many-body physics. They are the ground states of the isotropic 
Lipkin-Meshkov-Glick Hamiltonian [71] 



n 


£ (<*><*’ 


y y 


n 



i= 1 


(89) 


which describes n spins interacting through the two-body ferromagnetic coupling (A > 0) and h is the intensity of the magnetic 
field along the z direction. By tuning the value of h, the ground state of (89) varies along all \D%) for every k. Observe 
that the interactions in (89) are over all possible pairs, so it is worth mentioning that Dicke states of n qubits can also be 
well approximated for low n as ground states of a one-dimensional ferromagnetic spin-1/2 chain with nearest-neighbor XXZ 
interactions [72]. The states |D*) are genuinely multipartite entangled for any 0 < k < n. Recall that any //-qubit symmetric 
state (pure or mixed) is entangled if, and only if, it is genuinely multipartite entangled [73, 74]). These states are also generically 
genuinely multipartite nonlocal [75], Dicke states can be generated in experiments with current technologies (see, e.g.. Refs. 
[56, 76, 77]). 

In Ref. [52] we showed that the two-body Bell inequalities are capable of detecting nonlocality in some of the Dicke states, 
namely, \D^} with k = [N/2 J and k = \N/ 2]. Here, we show that nonlocality of any of the entangled Dicke states can be 
revealed with such Bell inequalities. In this section we present two classes of two-body Bell inequalities that are violated by any 
Dicke state | D%) for 0 < k < n. Note that \D°) = |0)® n and |D") = |1)®”, which cannot lead to nonlocal correlations, as they 
are clearly separable states. The first class, introduced in Theorem 10 detects all Dicke states | D*) for which 0 < k < 3n/10 
(or 7n/10 < k < n, by an appropriate symmetry of the Bell inequalities, consisting on renaming the outcomes of all the 
observables); the second one is introduced in Theorem 11 and it detects the rest, where 3n/10 < k < 7n/10. As it can be 
seen in Figure 5, the limits < have been taken very conservatively, and there is a clear overlap between the two classes, both 
detecting Dicke states with a number of excitations k around 3n/10 (or 7n/10). For the states | D%), the expectation value that 
a Dicke state takes with a symmetric two-body Bell inequality is straightforward to calculate via Theorem 4 and it is simply 

(D*\B\D*)=d k . 

Theorem 10. The inequality (28) defined with the following parameters 

/3 C = (1 + 2 k)[(n — 2k — l) 2 + n — 1], a = /3 = (1 + 2 k)(n — 1 — 2k), 7 = e = k, 6 = k + 1, (90) 

where 0 < k < (n — l)/2, is a valid Bell inequality for P®" for any n. 

Theorem 11. Let v := \ n/2\ — k. The inequality (28) defined with the following parameters 
• If n = 0 mod 2: 


P c = ' [n + 2(2i^ + l)], 


= 2 un(n — 1), = a/r 


7 = n(n — 1), 5 = n, e = —2. 


(91) 


• If n = 1 mod 2: 

/ 'n 


[n + 3 + 4 v(v + 1)], a = (1 + 2v)n(n — 1), (3 = a/n, 7 = n(n — 1), 6 = n, e = —2. 


(92) 


where 0 < k < [n/2\ — 1, is a valid Bell inequality for P® 


for any n. 
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A. Quantum violation of the Dicke states 

Here we analyze which are the optimal measurements to be performed on a Dicke state | D%). By Theorem 4, the computation 
of the expectation value (D^\B(ip, 9)\ D*) is immediate because it is precisely d k , as defined in Theorem 4. Note, however, that 
the state is fixed, so we cannot suppose that (D^\B(ip, 9)\D%) depends only on 9 — ip, since the Theorem 5 does not apply. Using 
the results stated in Theorem 12, one can easily generalize Theorem 4 to the use of any pair of measurements Mo = no - a. 
Mi = ni • <7. However, we have seen that (D^\B(no, fii)|£)„) depends only on the difference between the azimuthal angles 
that define n 0 , ft! in spherical coordinates, and that this expectation value is minimal when such difference is zero. That implies 
that there is no loss of generality in considering real observables, and the present form of Theorem 4 applies. For the case of 
inequality (90), i.e. with Dicke states with low number of excitations, we have that the expectation value (D^\B((p, 9)\D%) is 
a symmetric function of <p and 9. Numerically we observe that its minimum is obtained at ip = — 9 , a fact which significantly 
simplifies our analysis, and it allows us to obtain an analytically closed form for its maximal quantum violation for a fixed state 
| D%). After some algebra, by solving 


_d 

d 9 


d fc (-M) = 0 


we obtain the following non-trivial solutions: 


9* = ±2 arccos 



(93) 


(94) 


where Ci = k[n — (1 + 3fc)] and Ci = n 2 (2k + 1) — n(8k 2 + Ak + 1) + 2k 2 (1 + 4 k) and the two ± signs are independent. 
Using this optimal angle 9* we obtain the maximal quantum violation, given by 

AC 2 

d k (-9*,9*) = € 0{{n 2 ) 2 /n 3 ) = 0(n). (95) 

For the case of inequality (91) or (92), a proof that quantum violation exists for any n was already given in [52] for k = \n/2]. 
However, such violation is not the maximal one, as it was assumed that Mo = cr z in order to simplify the calculations. 


We have calculated numerically the maximal quantum violation of (90)-(92), and the results are reported in Figure 5 for 
n = 2 10 . 


VII. CONCLUSIONS 

In this paper we have constructed and classified Bell inequalities involving 1-body and 2-body correlators of many-body 
systems. We studied maximal quantum violations of these inequalities, as well as robustness of these violations and dependence 
with respect to various parameters involved. We paid particular attention to Bell inequalities that are violated by Dicke states 
that have physical relevance: on one hand they can be prepared in experiments, on the other hand they appear as ground states 
of physically relevant models, such as the Lipkin-Meshkov-Glick model. 

Let us conclude by listing several experimental setups in which nonlocality in many body systems may be tested using our 
inequalities: 

• Ultracold trapped atoms. Dicke states have been recently created in spinor Bose-Einstein condensates (BEC) of Rubid¬ 
ium F = 1 atoms, via the parametric process of spin changing collisions, in which two rriF = 0 collide to produce a 
nip = ±1 pair [56]. These recent experiments demonstrate the production of many thousands of neutral atoms entangled 
in their spin degrees of freedom. Dicke-like states can be produced in this way, with at least 28-particle genuine multi¬ 
party entanglement and a generalized squeezing parameter of 11.4(5) dB. Similarly, Rubidium atoms of pseudo-spin 1/2 
in BEC may be employed to generate scalable squeezed states (for the early theory proposal see [78], for experiments see 
[79]). Very recently non-squeezed (non-Gaussian) entangled states of many atoms [80] were generated in this way. The 
number of atoms used in these experiments are of order of thousands and larger. So far, these numbers and experimental 
errors and imperfections are too large, while the corresponding fidelities too small to detect many body nonlocality. In 
principle, however, it is possible to perform these experiments with mesoscopic atom numbers (say < 100, controlling the 
atom number to a single atom level (see Ref. [81] for the resonant fluorescence detection of Rb 8 ' atoms in a MOT, and 
Refs. [82, 83] for optically trapped spin 1/2 fermions). 
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FIG. 5. (Color online) On the left, the values of the maximal quantum violation, denoted \Q V \, of (90 - 91) for the Dicke state | D*) with 
n = 2 10 = 1024 qubits. The violation is normalized to the corresponding classical bound /3 C . The blue curve corresponds to (90) and the 
red curve corresponds to (91). The maximal quantum violation has been plotted in logarithmic scale (in dB) in order to better compare the 
two bounds; i.e. the actual plotted value is 101og 10 (—min v ,g dk/Pc), where d k is the quantum expecation value {D^\B(ip, 9)\D^)- The 
horizontal axis is cut at k = n/2 since the plot is symmetric with respect to k = n/2. The reason for that is that the same violation can be 
achieved by an inequality that results from renaming the outcomes of the measurements or, equivalently, by switching |0) <->• |1) for every 
qubit. It is clear from the plot that the two classes overlap, so these two classes cover all possible entangled Dicke states | D^). 

On the right, the optimal measurement angles tp* /it, 9* /-k for every k that lead to the violation shown on the left. For the class of inequalities 
(90) (j >* and 6* are plotted in blue and red, respectively. For the class of inequalities (91), the optimal measurement angles ij>* and 9* are plotted 
in yellow and green, respectively. Interestingly, although the plot suggests that ip* = n for the inequality (91), it is not the case. If it were 
exactly 7r, one could not achieve the optimal quantum violation. For this particular n, </>*|fc =280 ~ tt and rf>*\ k=511 « 1.00023-7T. Intuitively, 
the reason for this slight discrepancy is that d k can be expressed as d k = 2 n 4 cos(ip/2) 4 — 2 n 3 cos(ip/2) 2 (fe — cos# + (1 + 3fc) cosy)) + 
2n 2 k 2 cos 2 (ip/2)(l + 3cosy)) + o(n 2 k 2 ), and this expression can be made o(n 2 k 2 ) just by picking ip = 7r. However, one can do better 
by taking advantage of the fact that a small enough value of cos(y>/2) makes the whole expression effectively o(n 2 k 2 ) (this value of ip* will 
actually depend on n and fc), and the leading terms actually help in the minimization in this situation. 


• Ultracold trapped ions. Ultracold trapped ions with internal pseudo-spin ’’talk” to each other via phonon excitations, 
and under some conditions behave as spin chains with long range interaction. This was originally proposed in Ref. [84], 
using inhomogeneous magnetic fields, and in Ref. [85], employing appropriately designed laser-ion interactions. The 
pioneering experiments were reported in Refs. [86, 87]. While in the early theory studies [85, 88-90] spin interactions 
decaying with the third power of the distance were considered, it was experimentally demonstrated that management 
of phonon dispersion allows to achieve powers between 0.1 and 3 in the 2D arrays of traps [91], Recent state of art 
represents the work on experimental realization of a quantum integer-spin chain with controllable interactions [92]. We 
have studied trapped ion systems in relation to long range SU( 3) spin chains and quantum chaos [93], and trapped-ion 
quantum simulation of tunable-range Heisenberg chains [94]. In the latter paper we demonstrated that significant violation 
of the Bell inequalities, discussed in the present paper, is possible even for the ground states of the models with large, but 
finite interaction range. 

• Ultracold atoms in nanostructures. Yet another possibility concerns systems of ultracold atoms traps in a vicinity of 
tapered fibers and optical crystals (band gap materials). The experimental progress in coupling of ultracold atomic gases 
to nanophotonic waveguides, initiated by Refs. [95-97], is very rapid (cf. [98]). Early theoretical studies concerned 
remarkable optical properties of these systems (cf. [99-102]). Ideas and proposals concerning realization of long range 
spin models were developed more recently, and mainly in Refs, see [103-105]. 

• Cold and ultracold atomic ensembles. Last, but not least, one should consider cold and ultracold ensembles (for an 
excellent review see [106]), in which, by employing quantum Faraday effect, one can reach unprecendented degeees of 
squeezing of the total atomic spin (cf. [107, 108], and unprecendented degrees of precision of quantum magnetometry 
(cf. [109]). Note that in many concrete realisations the many body Bell inequalities derived in this paper require precise 
measurments of the total spin components, and their quantum fluctuations. Quantum Faraday effect, or in other words 
spin polarization spectroscopy, seems to be a perfect method to achieve this goal; note that in principle it allows also to 
reach spatial resolution, and/or to measure spatial Fourier components of the total spin [110], 
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Appendix A: Proofs 
1. Proof of Theorem 1 


Proof. Let us start from the if part and consider p = (a, 6, c, d) £ T„ \ dT n . This implies that all the coordinates of p are 
strictly positive, i.e., a, b,c,d> 1. We will show that p is not a vertex of P® " by explicitly finding a convex decomposition of p 
into vertices of P®". Let us consider a vector v = (1, — 1, — 1,1). Note that v is proportional to the fourth row of the Hadamard 
matrix in (29). Since the Hadamard matrix is orthogonal, all its rows are pairwise orthogonal. Consequently, for any A, p + Xv 
has the same value for n, So and Si, whereas Z(p + Xv) = Z(p) + 4A. Hence, the only coordinate that can change is Soi, where 
one obtains Sqi(p + Xv) = Soi(p) — 4A. Observe that for any Mi, M2- we have 

AH p{p + P2v) + p-2pip - Piv) = (ah + AH )<p(p)- (Al) 

If ah, M 2 > 0, then a convex decomposition of p{p) is given: 

V>(p) = —x— 4 - X—<m(p - W*') ( A2 ) 

Mi + M2 Mi + M2 

Notice that the choice ah = min{a, d} and ah = min{6, c} ensures that p + ah'!' £ 9T n and that p — p\v £ dT n . Since 
Mi, M2 > 0, this ensures that (A2) is a proper convex decomposition, so ip(p) Ext(P® n ). 

Conversely, for the only if part, let us consider p £ <9T„. If pip) were not a vertex of P® n , it could be decomposed as a 
convex combination of elements from Ext(P®"), namely {(/,;}(Lj with some k > 1. Since p is an invertible transformation in 
deterministic local strategies, for each qi there is a unique pi £ T ra such that p{pt) = qi- To see that p is indeed invertible in any 
local deterministic strategy, it suffices to observe that (So, Si, Soo,Soi, Sn) = (So, Si, (So) 2 — n, SoSi — Z, (Si) 2 — n), from 
which one can trivially obtain n, Si, So, Z and, from these, a, b, c, d by virtue of Equation (29). Denote then p, = (a^, bi, Ci, di)\ 
p(p) has then the following decomposition: 


v(p) = ^\v{Pi), 


(A3) 


i =1 


where 0 < Xi < 1 and the sum of the A* is 1. 

By looking at the third coordinate of Equation (A3), the fact that Soo(r) = (So(r)) 2 — n for all r £ T„ implies the following: 


(S 0 (p)) 2 - n = X A,[(S 0 (pi)) 2 - n) 


2=1 

2 


X x i s o{Pi) = X x i(So(Pi)) 2 ■ 


(A4) 


k 2=1 


2=0 


Treating Equation (A4) as a quadratic function in So(p m ), for some 1 < m < k we obtain 


^m(^m — l)(^o(Pm)) 2 + 2A m <So(p m ) ^2^0 (.Pi) + ( K$o(Pi) j — ^i(«So(Pi)) 2 — 0- (A5) 


i^m 




i^m 


Equation (A5) has a real solution So(p m ) £ R if, and only if, its discriminant is non-negative, which is equivalent to 

- 4A m X XiXj(So{pi) - S Q (pj)) 2 > 0. 
i<j- 


(A6) 
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Because all A,; > 0, (A6) is satisfied if, and only if, So(pi ) = So(Pj) for all h j 7 ^ m ■ However, this is true for every to, which 
allows us to conclude that it must necessarily be 6 >o(Pi) = <So(p) for all i. By applying exactly the same reasoning to Si, we 
finally obtain 

Sx(p-i) = S x (p), 0 < x < 1, 1 <i< k. (A7) 

The assumption that y>(p) ^ Ext(P®") leaves only the possibility that all <Soi(jPi) must be different. Since in any local deter¬ 
ministic strategy one has Soi(r) = So(r)Si(r) — Z(r), Equation (A7) implies that all the Z(pf) have to be different. The fourth 
coordinate or (A3) then reads 


k 

a — b — c + d= Ai(<Zi — bi — Ci + dj). (A 8 ) 

i —1 

Equations (A 8 , A7) together with the fact that a + b + c + d = ai + bi + + di = n imply (by appropriately adding and 

subtracting them) the following: 


k k k k 

CL ^ b ^ C ^ ^ d ^ ^ (A9) 

z— 1 z=l z=l i= 1 

which is equivalent to p = 2 ^i=i ^iPi\ i- e - P is a convex combination of elements of T n with the same weights as the convex 
combination of <p(p). Finally, if p £ <9T„, then one of its coordinates must be 0. However, if any of its coordinates is 0, 
Equations (A7, A 8 ) and a + b + c + d = ai+bi + Ci+di imply that p = Pi Vi, contradicting the fact that (A2) is a proper 
convex decomposition. This contradiction comes from the fact that we assumed that ip(p) was not a vertex of P ®". Hence, 

p(p) e Ext(P®"). □ 


2. Proof of Theorem 2 


Proof. We consider the function / defined as 


^ 2 2 
I ■= 2 [ n ( x + y ) 2 + ± x ) 2 - 1] + x[ay ± (x + y)]5 0 + yySi + 2 —S 00 + axySoi + y <Sn 

and we wish to show that min TP a 1 = 0. For all local deterministic strategies, I has the form 

2 2 2 
/ = — \n(x + y) 2 + (ay ± x) 2 - 1] + x[ay ± (x + y)]S 0 + yySi + y (So - n) + axy (S 0 Si - Z) + y (S? 

= - ( x 2 Sq + 2crxySoSi + y 2 S 2 ) + xyn — axyZ + -(cr/i ± x) 2 + (ay ± :r):rSo ± xySo + yyS\ — - 

= - (xSo + aySi) 2 + - [(cry ± x) 2 + 2(ay ± x)xSq + 2(ay ± x)aySi] =F axySi ± xySo — axyZ + xyn - 

1 2 1 
= - (xS o + crySi + ay ± x) + xy (±S 0 T crSi - aZ + n) - 

The condition I > 0 can be recast now as 


(A10) 



1 

2 


(xiSo + ay Si +ay± x) 2 + 8xyr > 1, (All) 

where r := (±<5>o T crSi — aZ + n)/4. We note that r results into one of the variables a, b, c, d which depends on the choice of 
the signs a, ± we make. Indeed, one has the following values for r: 


cr\± 

+ - 

+ 

b c 

- 

a d 


(A12) 


The left-hand side of (Al 1), denoted I, is always non-negative, because x, y £ IN and r > 0. We wish to show however that 0 
cannot be achieved, thus proving that I is indeed tangent to P®". Note that if r > 0 or xSq + ayS\ + ay ± x 0, inequality 
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(A11) is trivially satisfied. Thus, the set of points in which 1 = 0 could be possible needs to be a subset of Jo, where Iq is defined 
as 

Io := {(a, b, c, d) £ T„ : r = 0, xSq + ayS\ + ay ± x = 0}. (A13) 

Observe that when we take (a, b, c, d) to be continuous, then Iq is a line lying on a facet of dT n . We shall now see how the 
conditions of Theorem 2 ensure that Iq = 0. Let us discuss the parity of xSq + ayS i + ay ± x. Note first that for all m £ 7L, 
m = m 2 = —to mod 2. Then, it follows that Sq = Si = n mod 2. If n is even, xSq + aySi + ay±x=y + x mod 2. 
But y + x = 1 mod 2 by hypothesis. Otherwise, if n is odd, then xSq + crySi + ay ± x = x + y + y + x = y + y mod 2 
and in this case y + y = 1 mod 2 again by hypothesis. Hence, Iq = 0, which implies (All), which in turn proves that (49, 50) 
is a valid Bell inequality, tangent to P®" for any n. 

□ 


3. Proof of Theorem 3 

Proof. We will prove the theorem by explicitly solving the Diophantine equation 

-\n{x + y) 2 + (ay ± a;) 2 - 1] + x[ay ± (x + y)]S Q + yySi + ^Sq 0 + axyS 0 i + ySn = 0 (A14) 

over the integers. As discussed in the proof of Theorem 2, (A14) is saturated on r = 0 and ( xSq + ayS\ + ay ± x) 2 = 1. 
Let us rewrite the last condition as xSq + ayS\ + ay ± x + r = 0, where r £ {—1,1}. In the notation introduced in (55) 
this condition reads K r (r) ± 2 (xs — yt) = 0, where K r (r) is defined as K r (r) := ±n(y — x) + ay ± x + r. Let us notice 
that the assumptions of Theorem 2 guarantee that K r (r) is even: Indeed, if n = 0 mod 2, then K r (r) =y+x+l=0 
mod 2 (because the assumptions on Theorem 2 guarantee that y has opposite parity to x when n is even); otherwise, if n = 1 
mod 2 then K r (r) = y + y + 1 = 0 mod 2 (because y has opposite parity to y when n is odd). Hence, we can define 
K r (r) \= K r (r)/2 and we just showed that K r (r) £ Z. 

Now, the condition K r (r) ± 2(a;s — yt) = 0 is equivalent to K r (r) ± (xs — yt) = 0, which we can solve for t: yt = 
±K r (r) + xs. This equation, when taken modulo x , reads yt = ±K r (r) mod x. Since gcd(a :,y) = 1 by hypothesis, this 
equation has a solution t = ±y~ 1 K r (r) mod x , where y~ l is the inverse of y in the group of integers modulo x, typically 
denoted Z., ; , (such an inverse exists if, and only if, gcd (x,y) = 1). In practice, computation of y~ l is done via the Bezout 
identity, which says that for any pair of integers x, y £ Z there exist (not unique) p, q £ Z such that px + qy = gcd(a;, y). 
In the case that gcd(a:, y) = 1 when taking Bezout’s identity modulo x we obtain qy = 1 mod x\ because of the assumption 
gcd(x, y) = 1 of the Theorem, the inverse of y modulo x is well defined and we can write q = y~ l mod x. 

Let us now solve Equation (A14) for the other variables s, u, but we will keep the following form for f, in order to count the 
total number of solutions: t(r) = to(r) + kx , where k £ Z and 0 < t 0 (x) < x. Then, from K r (r) ± (xs — yt) = Owe solve 
for s and obtain s(t) = sq(t) + ky, where so(t) := (=F K r (r) + ytQ(r))/x. Note that so( T ) is well defined integer number, a 
fact that can be shown by directly solving xs(t) = (=F K r (r) + yto(r)) + kxy , because yt(r) = yto(r) = ±K r (r) mod x. 
Now, the last variable u is directly obtained from the normalization condition n = r + s + t + u. Then, writing u = n — s — t as 
u(t) = u 0 (t) — k(x + y), where u 0 (t) := (nx + K r (r) — (x + y)to(r))/x, which is integer by the same argument as s 0 (r). 
To sum up, the family of points (r, s, t , u) £ Z 4 for which (A14) is fulfilled is 

IJ {[0,s o (r),to(T),uo(T)] + k[0,x,y,-(x + y)] : k £ Z} (A15) 

TG{ —1,1} 

Geometrically, this corresponds to alternating points in a zig-zag pattern along two parallel lines which lie on the facet r = 0 of 

T 

R n- 

Let us now study which is the number of solutions that really belong to T n ; i.e. those for which r, s.t. u > 0. Since we have 
the family of solutions indexed by k, we just have to count how many fc’s are available. The condition / > 0 leads to to + kx > 0, 
which for k £ Z means k > [— t o /a?]. Now, taking into account that f 0 is chosen to be 0 < to < x, which is equivalent to 
0 > —to/x > —1; equivalently [— to/x) = 0. Thus, k > 0. The condition s > 0 becomes sq + ky > 0, which for k £ Z is 
k > f—so/2/1 ■ Finally, the condition u > 0 is uq — k(x + y) > 0, which for k £ Z is k < \uq/(x + y) J. 

Hence, for each r £ { — 1,1}, the number of solutions belonging to T n is given by |{fc £ Z : k > max{0, [— so(r)/y ]}, k < 
\uq(t)/(x + y)\ }|, which is max{0, \ uq(t)/(x + y) J — max{0, {— So(r)/y]} + 1}. Note that we put the first maximum just to 
ensure that there is always a non-negative number of solutions, denoted Ng. Finally, we obtain Ng by summing this expression 
over the possible values for r: 

Ng = ^ max < 0, 

r=± l *- 


uo(t) 


— max < 0, 


-So(t) 


+ 1 


(A16) 








28 


□ 


4. Proof of Theorem 4 


Proof. Recall that the Bell Operator can be decomposed into a sum as 


B(<p, 0) — a <So + PSi + ^Sqo + 5S 0 i + fSii + /3 c l2~ (A17) 

= AS Z + A' S x + BS ZZ + CS XX + DS XZ + f3 c I 2 ™, 

where and S fcl = E^ with k, l £ {at, 2 }. By projecting everything into the space PLj of (59) we 

obtain the analytic expression of every block Bj(tp, 9), as the projection of ,S7 and .S'/,. / is given by Theorem 12. Simply by 
joining terms we arrive at the form (66). □ 


5. Proof of Theorem 5 

Proof. By a direct check one finds that for any c £ R, U(c)Aio{ip)U(cy = A4o(ip + c), where (7(c) is given by Eq. (75). 
Obviously, the same holds for A4i (6). Consequently, 

B(<p + c,9 + c) = [U(c)f n B(ip,6)[U(cY}® n , (A18) 

which in particular means that the spectrum of B(ip, 9) is invariant under mutual translations of both its arguments, and therefore 
it depends only on the difference p — 9. □ 


6. Proof of Theorem 6 

Proof. We have to make C = 0 with C being defined in (73). Hence, 9± is the solution to C = 7/2 sin 2 <p + <5 sin sin 0 + 
e/2 sin 2 9 = 0. This means that 


sin ip = 


—26 sin 9 ± \J 4i5 2 sin 2 9 — 47 c sin 2 9 


27 


-6 ± 


7 e 


sin0. 


—5 sin0 ± \/6 2 - 7e| sin0| 

7 


(A19) 


Note that in the third equality the absolute value is omitted due to the freedom of choice in the sign ±, which can be absorbed 
into the sin 9. 

Let us denote 


z~ 


-5±VA 

7 


A := 6 2 — 7e. 


Then, we wish to solve for 9 the equation sin(0 — ft) = £ sin(0) for k £ R. We can assume that n Y miT 1 m £ Z, otherwise 
the equation is trivially satisfied by setting 0 = 0. But then A4o = ±A4i, which is a degenerate case. Then, by expanding the 
sin(0 — ft), one obtains 


£ 


sin 9 cos k — cos 9 sin k 
sin 9 


cos k — 


sin k 
tan 9 ’ 


which leads to 


9 = arctan 


sm k 

COS ft — £ 


= arctan 


7 sm ft 


7 cos ft 


■ 6 zb v'A 


□ 
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7. Proof of Theorem 7 


Proof. Let us first recall a known result from probability theory: The non-central moments of the gaussian distribution are given 
by ' 


£„,„(**) = J 


°0 e ~ (x—fj,) 2 /2(7 

X -; - - 


dx = a k/2 H k (^a 1/2 ) (k = 0,1,...), 


where H k {x) are the Hermite polynomials, which can be defined as 


H k (x) = e 


- p -* 2 / 2 _ 


dx k 


p * 2 /2 


So, we have 

{E MitT (x fe )} fc >o = {1, n, /i 2 + <7, p. 3 + 3/rcr, /t 4 + 6 ^ 2 ct + 3cr 2 , /r 5 + 10^ 3 cr + 15/nr 2 ,...}. 
Additionally, we shall make use of the two following facts in the proof: 

• Let ci, C 2 £ R. Then 


(A20) 


roo g — (x—ii+c 1 ) 2 /4cr — (a: — /i-)-c 2 ) 2 /4<r 

X fc -- 

—oo v27T(J 


da; = e- (ci - C 2 ) 2 / 8 ' T E A 1 ,, (T (a: fe ), 


(A21) 


where jj! = /x — (ci + cf)j2. This can easily be verified by using the following expression 


(x — a) 2 + {x — b) 2 


x — 


a + b 


+ 


a — b 


with a = n — ci and b = // — C 2 - 
• Let us consider the family of states 


1 n 1 n 


o-(fc-M) 2 /4cr 


k=0 


Nn 


\D k n ), 


k =0 


where /u = n/2 + A/{2B — C) and e 1 /27r <C ct <C n. Then, one finds that limn-j.oo A/" n = 1 because for large n the 
following holds: 

1 " e ~(k~ l i) 2 /2a i rn - (fc- M ) 2 /2cr 1 


fc =0 


AC 2 ' 


Observe that the choice of er cannot be too small: otherwise the first approximation (the sum by an integral) would be 
invalid; o cannot be too large either as the second approximation (integration over R instead of the interval [0, n]) would 
be incorrect. This is the reason why we added the constraints e _1 /27r <C a -C n. 

Let us then assume that n is large enough so that we can omit the factor Af n for practical purposes. To calculate (t/v, \B n / 2 {}p{0), 9)\f) n ), 
let us first express d k and u k conveniently as polynomials in terms of k (note that C = 0): 

• d k = 2 Bk 2 — 2(A + nB)k + \fl c + n{A — B/2) + n 2 B/ 2], for 0 < k < n. 

• u k = Uk'+m-i = Uk' = {A'—2k'D)myJ 1 — for —m+1 < k' < to— 1, where m = (n+l)/2 and k' = k— (n— 1)/2. 

By means of the Taylor expansion of \/l x' 2 around x - 0, 


°° 111 

\/\ — X 2 = Y C\X 21 = 1- x 2 - x 4 - x 6 — 

v 9 8 16 


1=0 


128 


+ o{x 1U ), 


we can express v! k as a polynomial in k'\ this is a good approximation as long as |fc'| <C to: 

12 k' 4 k' 6 5k' 8 


u k > = ( A' - 2k'D) ^ - 


8 to 3 16?n 5 128?n 7 
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• Vk = 0 for 0 < k < n — 2. 

Then, the expectation value of the elements of the diagonal of B n /i{<p(Q), 9) is given by 

r n e -(k-ii) 2 /2a 


pn 

/ dk- 

k =0 do 


7TCT 


r°° p-(fc-A*) / 20 - 

-dfc~ / d k - - -—dfc 

J —oo V 27TCT 


= 2BE^ a (x 2 ) - 2(A + nB)E^(x) + [/3 C + n(A - B/2) + n 2 B/ 2]E #1 , <r (l) 
= 2H(/x 2 + cr) - 2(A + nB)p + \/3 c + n{A - B/2) + n 2 B/ 2] 


& S 


= 0 ■ n + [ - -— n + 2,Bo — ——, 


n 2 


A 2 


2 B' 


(A22) 


where we have applied (A20). 

The contribution to the expectation value for the elements that correspond to the off-diagonal terms of B n /2(<p{9), $) is given 
by (let us write in short c = A/2B) 


n—1 m—1 

i,(")„/,(” 


X'n" ) n+i u k= X 

k'=— (m— 1) 

pm— 1 g— (fc'+m— 1— fi) 2 /4a— (k 1 -\-m—fi) 2 /4a 


U k ' 


' — m+1 


dfc' 


7T(T 


p- 1 e -(fc'-(c-l/2)) 2 /4<T-(fc'-(c+l/2)) 2 /4<T 

/ ™ k ' -771=- 

/ — m+1 V27T(J 

,oo g—(k' — (c—l/2)) 2 /4<7—(k'—(c+l/2)) 2 /4a 

/ -771=- 

/ — oo y 27T(J 


dfc' 


dfc' 


,/2 k ,4 


/ k' 2 


k' 6 5k' 8 


8 m 3 16m 5 128 itA 

e _( fe '_ (c _l/2)) 2 /4<T-(fc'-(c+l/2)) 2 /4a 


X 


dfc' 


/ oo / 7/2 

.+ - 2k'D) ( m — 2m — 


,/2 i ,/4 


V Z7T<T 

fc ' 6 5fc ' 8 


8 to 3 16m 5 128m 7 


v/ 27rcr 


-dfc' 


= e 1/8ff (A'mE clT (l) - 2£>mE c ff (x) - ^-E C(7 (a; 2 ) + —E CjCr (x 3 ) 

\ zm m 

,(A! - 2Dc A'~2Dc (c 2 + a)A'- 2Dc(c 2 + 3a) 2 ,.\ 

= e- 1 '*" (—j—» + —j— - ! - i —+-STl + 0<a " >) ■ 

where we have applied (A20 and A21). Hence, the expectation value (+|B 7l / 2 (+0), 9 )|+) is given by 

n 2 n—1 

(4>n\B%(ip(9), 601+) = X (^X) dk + 2 X 




n 2 


( A ' - +))" + (+ - 4 + e ~ m " l' 4 ' - +)) 


□ 


8. Proof of Theorem 8 

Proof. Since the state acts on the symmetric space, without loss of generality we will compute its partial trace after forgetting 
about the last n — d subsystems. In general, when p acts on the Hilbert Space of n qubits (C 2 )®", its parcial trace after forgetting 
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about the last n — d subsystems is given by 


(Tr n _ d (p)) 


lQ...ld-l 


E 

0 <l d ,...,ln 


Iq , • •. i^d— 1 > • • • i^n— 1 


Now, let us expand ps to the full Hilbert Space: 

(p) j = (ppsp t ) ] = Y Pk(ps)Wi = E \ ps ^ 1 s ( k ~ 1*1)^ “ IjI)> 

J J J ^— J / (n\ (n\ 

0 <k,l<n 0<k,l<n \l \ k ) \i) 

where p is the same matrix for the change of basis as the one defined in the proof of Theorem 12, but in this case m = 0. 
Joining the last two expressions and introducing the notation i ' := id ■ ■ ■ i n - 1 and j' := jd ■ ■ ■ jn-i, we have 


(Tr„_ d (p))! ; = (Tr n-d(ppsP T ))\, = Y Y 


(Ps)f 


= E E 


T 0<fe,i<n ;/ (l) (7 
( Ps)f 

0<k,l<n V (£)(?) 




1 




= £ -&S=<s(fc - |i'| - 1 + |j'l) E S ( k ~ I 1 '! “ I 1 '!)’ 

0<M<n t/m F 


(A23) 


where in the third equality we have used that Kronecker deltas fulfill S(a — b)5(a — c) = S(a — b)5(b — c). The expression 
X^irc)(fc — i 7 — |i'|) is equal to (fe"L|^|)^[|i'|,|i'|+n-d] (A:), where I[ a ^\{k) is the indicator function: it outputs 1 if a < k < b and 
0 otherwise. The reason behind that equality is that we are counting how many indices i with the first d bits already set and with 
weight |i' | have total weight k. Hence, we have to choose from the remaining n — d bits, indexed by i 7 , that they have weight 
k — |i 7 |. Of course, this is possible only if |i 7 | is not already greater than k and it is not smaller than n — d k\ that is the reason 
why the indicator function appears, so that |i 7 | < k < |i 7 | + n — d. We can then write 


(Tr„_ d (p)) 


j' 


E 

0</e,Z<n 


(ps)? 

vEk?) 


S(k- |i'| 



HjWI n — d^ 

v'oo-i.-Wi) 


-f[|i'I,|i'|+n—d] {k) 
-A[|i'|,|i'|+ri— d]{k) 


E 

\i' \ <.k<.\i' \ -\-n—d 


(As)fc— |i'|+|j'| f n-d^ 
11,17 


Finally we apply the change of variables m := k — |i 7 | in order to obtain 


ln—d\/ \m+|i'| 

(Tr„_ d (p))i: = Y, / ' ^ 


0<m<n—d 






which completes the proof. Observe that the expression only depends on the weights of the indices, so it is permutationally 
invariant, as expected. □ 


9. Proof of Theorem 9 

Proof. We shall make use of the formula provided in Theorem 8 and use the same kind of approximations as in the proof of 
Theorem 7. We shall use the definition for the function / given in Eq. (83) and its relation (81). We shall also adopt the notation 
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(A20). The simplest elements of p 2 to begin with are its diagonal elements: 


(P2)oo - fo 2 /0> x > 2 , °> 0 )V>x da; = n(n- 1) So 2 (n - x)(n - x - l) e da: 

- n(n-l) ( n ( n “ ^ ( 2n - l)R^,*(x) + E^ CT (x 2 )) , 

/(n,x,2,l,l)^ +1 dx = Jo~ 2 (n-x~l)(x + l) 


(P2)ol - fr 


(P2)\\ 


n—2 
0 


/2<t 


— n(n-i) (( n + ( n 2)E m _i j(J (x) E /i _i jtT (a; 2 )), 

~ r 

— Jo 


y/27T(T 


dx 


f (n, x, 2, 2, 2) , 0 2 +2 dx = /” ^(x + l)(x + 2) e U ^ ->1 '''" dx 


— n(n-l) ( 2 ®>-2 ,<t( 1) + 3E M _ 2) o-(x) + E M _ 2l(7 (x 2 )) . 


y/2TT(T 


Note that (^ 2)01 = (^ 2)10 = (P2)oi = (^ 2 ) 10 * Observe as well that Trp 2 does not depend neither on /1 nor a and it is 1. Now 


we move 


to the off-diagonal elements, starting with (/? 2 )ii> which is equal to (^ 2 ) 


11 

oo- 


fo)?? 


- So 2 0 , 2)ip x il; x+ 2dx 

= nfn^l) IS 2 X )( n - X ~ 0( X + + 2 ) ' 


(x-m ) 2 -(x-/x+ 2 )^ 


n(n—1) JO 

- ^1/2)(x + 3/2)^ 

= 1 / 2 -»)(» + 1 / 2 )* 


(x-( M -l)) z 
_ 2 <T 

y/2'Kcr 
-(y-v) 2 / 2 cr 


\Z2na 


-dx 


-dx 


y/2'Kcr 


dy 


n(n— 1) 


((^+ 1 ) Emi(T ( 1 ) + nE Mi<r (aO - E^x 2 )) . 


Here we have used the approximations y/(n — x)(n — x — 1) ~ (n — x—1/2) and sj{x + l)(x + 2) ~ a:+ 3/2 and the change 
of variables y = x + 1. Finally, we consider the cases (p 2 )oi = (^ 2)10 = (^ 2)00 = Moo and (^ 2)11 = (^ 2)11 = (^ 2 ) 0 ! = 
G° 2 )io- To this end, let us recall that y = n /2 + A/2B. We shall denote c = A/2B and we shall also define m = [n + l)/2 
and /a' = y — 1/2. 


(+ 2 ) 


_ 9 _ 9 _ (x-^) 2 -(x-^+1) 2 

-fo f(n, x, 2,0,1) AA+idx = ^ny/ 0 " (n - x - 1)^^ - x)(x + l) 5 - 4<T 4<T -dx 

S-^n-x- l)y/(n — x)(x + l) e ^ dx = /^(m - 1 - y)my /1 - fa/m) 2 d V 


n(n—1) 

ff+t) SZJ 171 - 1 - y) m EZo c i(y/ m ) 


21 e 


-(y — c) z /2a 


e -!/8 

n(n 
p - 1 / 8 


y/2/KG 


dy 


f-oJ™ - 1 - y) TO ( X - i - 8 +* - 0 ((t//m) 6 )) 


,-(w- C ) / 2 <r 
y/2ncr 


dy 


- ( TO ( TO - 1 ) e c,<t(1) - toE c , ct (x) - ^E c , CT (x 2 ) + 2 ^E CjCT (x 3 ) - ^/E CjCT (x 4 ) + 8^E C!ff (x 5 )) 


e ~ 1 /SfT 

n(n— 1) 

e" 1 / 8 " 

n(n— 1) 


(n+l)(n-l) 
4 


++ 1 ) 


c — 


2(n+l) 


(c 2 + er) + ^(c 3 + 3ccr) + • • • ) 


In 2 - 9-71 - 2c z +2c+l+2a 
4 n 2 1 4 


C +(7 + C J +3C(7 


+ • 


)• 


where q are the coefficients of the Taylor expansion of s/1 — x 2 at x = 0, we have used the change of variables y = x— (n— 1)/2 


and the fact that 


1 


n — a 


1 

n 


E(«/ n ) fc 

fc =0 


(A24) 


in order to obtain an expression purely in powers of n. 


(P 2) 01 


11 


- /o" 2 f(n, x, 2,1, 2)t/j x+1 ^ x+2 dx = „ (n 1 _ n / 0 " 2 (x + l)v / («- ® - l)(x + 2)^ 

/s „ ,_ (»-(+-i)) 2 

3T) IZo( x + 0\/( n ~ x - 0( x + 2)--7^-dx 


4cr 4cr 


\Z2tvg 


dx 


e — l/8a 

n(n— 1) 
-1/8 


y/2-KG 

S SZo( m - 1 + y^V 1 - (y/ m ) 2 e ~ fv JS^ /2a dy 

TTT f-oo( m - 1 + 2/) w ( X - 2 ^ ~ sFH ~ 0((y/m) 6 


n(n— 1) 

p-l/S^x 


0 -{y-c) z /2a 


dy 


2m 2 8m 4 ^\\y/"“> JJ v ^r 

(m(m - l)E Ci<T (l) + mE v (x) - n ^E C:a {x 2 ) - 2^E c , CT (x 3 ) - f^E CyCr (x 4 ) - gi3E c , CT (x 5 )) 


i(n—1) 


( (n+1) l ra 1} + - 2f++I) (c 2 + °) - ^+t(c 3 + 3cct) + • • • ) 

g-i/So- / lm2 ( 2 c 2 +2c+1+2o- , c 2 +<t-(c 3 +3c<t) , / 

n(n-l) ^4 n 2 n 4 ^ n 4 / ’ 
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where we have used the change of variables y = x — (n — 3)/2. 
Joining all terms, we obtain that p 2 can be written as 



( 

/ l e-l/8 a 

e-i / 8ct e -1 / 2<T \ 


f —(2c + 1) 

—ce -1 / 8cr 

—ce 1 / 8cr 

e -1 / 2cr \ 

\ 

1 


e -t/8<r X 

1 e -1 / 8<T 

n 2 

—ce -1/,8<T 

0 

0 

ce” 1 / 8 " 7 

n / s 

P2 — , n x 

n\n — 1) 


e -t/8<r X 

1 e -1 / 8<T 

T + 

-ce- 1 /^ 

0 

0 

ce” 1 / 8 " 7 

2 + °( n ) 


\ 

y e -1 / 2<T e- x ' % ° 

g-t/ScT ]_ J 


\ e -1 / 2<T 

ce" 1 /^ 

ce -1 / 8<T 

-1 + 2 c) 

/ 


Since a £ Oin 1 / 2 ), we can neglect the exponential terms and we obtain the form given in Eq. (84). 


□ 


10. Proof of Theorem 10 


Proof. We begin by proving that the choice of coefficients (90) defines a valid Bell inequality; i.e ., (28) is non-negative on all 
points P®". To this aim, we shall express the inequality (28) in terms of the variables a, b, c, d and finding the minimum to its 
left hand side, while imposing the conditions a , b,c,d> 0. If we denote by F(a, b, c, d) the left hand side of (28), then we find 
that it does not depend on the variable d, and, as a quadratic function of a, b, c, it can be written in the following form: 



( 4fc(l + k)(l + 2k) 

2(1 + 2 k) 2 

—k( 3 + 4 k) 

—fc( 3 + 4 k) 

°\ 


/ 1 \ 


2(1 + 2k) 2 

4(1 + 2 k) 

2(1 + 2k) 

2(1 + 2 k) 

0 


a 

(1 ,a,b, c,d) 

—k{ 3 + 4 k) 

2(1 + 2 k) 

2k 

2(1 + k) 

0 


b 


—k( 3 + 4 k) 

2(1 + 2 k) 

2(1 + k) 

2k 

0 


c 


\ o 

0 

0 

0 

o / 


\d) 


We aim to prove that F > 0 on all vertices. The condition for extrema is d a F = d b F = d c F = 0, which reads 

d a F = 4(1 + 2k) [2a + b + c - (1 + 2k)) = 0, 

d b F = 2[2 b - (4 k + 3)]fc + 4c(l + k) + 4a(l + 2k) = 0, 

d c F = 2[2c - (4 k + 3 )]k + 46(1 + k) + 4a(l + 2k) = 0. 


(A26) 


(A27) 

(A28) 

(A29) 


We also have to take into account the boundary, so there are 8 cases to consider, depending on the case if the inequalities 
a. b.c > 0 are saturated or not, which are summarized in the following table: 


II 

o 

Cr 

II 

O 

o 

II 

o 

( a,b,c) 

F(a,b, c) 

(Y,Y,Y) 

(Y,Y,N) 

(Y,N,Y) 

(Y,N,N) 

(■ N,Y,Y) 

( N,Y,N) 

(■ n,n,y) 

(N, N, N) 

(0,0,0) 

(0,0,3/2+ 2 k) 

(0,3/2 + 2k, 0) 

/pv fc(3+4fc) /c(3+4fc) \ 

W 2(1+2 k) ’ 2(1+2 k) > 

(1/2 + k, 0,0) 

(k/2, 0,1 + fc) 

(fc/2,1 + fc, 0) 

The system is incompatible 

4fc(l + fc)(l + 2k) 
-k/2 
-k/2 

fc(8fc 2 +llfc+4) 

l+2fc 

— (1 + 2k) 
k 2 
k 2 

N/A 


(A30) 


Thus, in order to prove F(a, b, c) > 0, one only needs to look at the cases where F(u. b, c) is negative. Since k > 0, it is 
immediate to see from (A30) that this occurs when exactly two of the variables a, b, c are equal to zero. Since in these cases 
F is just a quadratic function of the remaining nonzero variable, we only need to check two cases, which are the neighboring 
integer values to the optimal (real number) given in (A30) for this nonzero variable. When exploring the case for which this 
nonzero variable is rounded below (cf. the case a > 0 in table (A31)) we must also take into account the case when the value of 
the variables set to zero becomes 1 or we might miss some solutions. The following table summarizes the strategies for which 
F = 0 (the rest are F > 0 so they are not listed) thus proving inequality (90). More precisely, from c > 0 one finds the first two 
rows of (A31), from b > 0 one obtains the third and fourth row of (A31) and the remaining ones are obtained from a > 0. 

Observe that the fact that (90) is symmetric under the exchange of measurements translates in the following: if a vertex 
exists for which F = 0, then by swapping b and c (recall that these variables corresponded to (A4<). ) = (+,—) and 

(A4o, A4i) = (— , +), respectively) we obtain another vertex for which F = 0. Equivalently, by swapping the values of So and 
Si, the same happens. 
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a 

b 

C 

d 

So 

Sr 

Z 

n 

0 

0 

1 + 2 k 

n — (2 k + 1 ) 

—n 

2 ( 2 k + 1 ) - n 

n — 2(2 k + 1 ) 

n 

0 

0 

2 + 2 k 

n — (2k + 2 ) 

—n 

4 (k + 1) — n 

n — 4 (k + 1) 

n 

0 

1 + 2 k 

0 

n — (2k + 1 ) 

2(2k + 1 ) — n 

—n 

n — 2(2k + 1 ) 

n 

0 

2 + 2 k 

0 

n — (2k + 2 ) 

4 (k + 1) — n 

—n 

n — 4 (k + 1) 

n 

k 

0 

0 

n — (k + 1 ) + 1 

2k — n 

2k — n 

n 

n 

k 

0 

1 

n — (k + 1 ) 

2k — n 

2 (k + 1 ) — n 

n — 2 

n 

k 

1 

0 

n — (k + 1 ) 

2 (k + 1 ) — n 

2k - n 

n — 2 

n 

k + 1 

0 

0 

n — (k + 1 ) 

2 (k + 1 ) — n 

2 (k + 1 ) — n 

n 

n 


(A31) 


□ 


11. Proof of Theorem 11 

Proof. We begin by making a general observation. For any vertex in P®", the following inequalities hold: 

- n < S 00 ,Su < n(n - 1), |<Soi| < n(n - 1), |<S 0 |, |<Si| < n. (A32) 

Hence, the dominating term in (91 — 92) is the one corresponding to 7<S 0 o/2. which is 0(n 4 ), while the rest are 0(n 3 ). Hence, 

in order to minimize the value of the Bell inequality, a necessary condition is that the term containing <S>oo should be small, at 

least not bigger than 0(n 3 ). As <Soo = (So ) 2 ~ n in the vertices of P®", this suggests to use So as a parameter and minimize 
the left hand side of (28) with the coefficients given in (91 — 92) for a fixed value of Sq. This leaves us with just two variables 
to consider in the optimization, because n and So are fixed. We pick, for example, b and d, so that we have 

_ ri + So _ ,,, n — So 

a — a(b) = —--o, c— c(d) = —-- d. (A33) 

Equation (A33) enables us to rewrite Si and Z as 

Si = n-2(6 + d), Z = S 0 — 2(6 — d). (A34) 

And the constraints a > 0, c > 0 impose bounds on b and d. which together with their non-negativity become 

0 < 6 < n ± 5 °, 0 <d< n ^ S °. (A35) 

~ ~ 2 2 

We shall proceed in the same spirit as in Theorem 10, and we start by focusing on the cases for even n and odd n, which have 
a similar proof: 

• Even 7i. 

We start by considering the Bell inequality (91) as a function of b and d, with n and So treated as parameters. We denote 

tt In,So d). 


- Case a = 0. Note that it is equivalent to b = (n + So)/2. Then we have 

In,So ( n+ 2 S ° 'dj = -4d 2 + h{y,Sa,n)d + f 0 (v,S 0 ,n), (A36) 

where f are expressions that do not depend on d. Equation (A36) is quadratic in d, and it fulfills that its second 
derivative is negative for all d, so Eq. (A36) has only one maximum. Therefore, its minimal value is attained either 
at d = 0 or at d = (n — Sf)/2. The first case does not lead to the optimal solution (See [52] for the case v = 0, 
where the same argument follows), so we shall omit it here. The second case produces the following result: 

In,S 0 = Q + 2 ^ + 2 ^ 2), 


(A37) 
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which is minimal for So = 1 — 2v. However, since we are in the case where n is even. So must be even as well, 
which now is not the case. So we look for the closest even integers, which are 2 — 2v and —2v. It is immediate to 
see from (A37) that the value of I is precisely 0 on these points. Hence, we find the following values of (a, 6, c, d) 
that saturate I: 

\ 0 ’2~ P,0, 2 +P ) , (0, 2 + g ~ 1 + v ) ■ (A38) 

- Case b = 0. In this case we need to consider 

4,s 0 (M) = + 9 i(y,So,n)d + g 0 {v,s 0 ,n), (A39) 

for some expressions g, that do not depend on d. Again, we see from (A39) that / has only one maximum, so its 
minimum must be either at d = 0 or at d = (n — 6>o)/2. If d = 0, we have that 

In, So (0,0)= Q (So + 2 v) (So i 2 n A 2 ), (A40) 

which is minimal for So = — 1 — 2v. Again, as <S 0 must be even, we consider two cases. So = —2 — 2v and 
So = 2v, for which Equation (A40) immediately guarantees that I = 0 on these points. We then obtain the 
following values of (a, b , c, d): 

/ n n \ fn n \ 

— 1 — 0, — + 1 + i/, Oj , ~ v i 0, — + v, Oj . (A41) 

If d = (n — So)/2 then we obtain 

In, So (o, = (n - l)(So + 2 v)((n + 2 )S 0 + 2nv)/2, (A42) 

which is minimal for So = —2which is not an integer. The closest even integers to So are —2v and —2v + 2 
and in this case, only —2v leads to I = 0 , which corresponds to the following coordinates (a, b , c, d): 

Q-t/,0,0 (A43) 


The points (A38, A41 and A43) comprise all the vertices of in which I is minimum, with corresponding value 0, 
showing that (91) is a valid Bell inequality, as there are no more cases. 

• Odd n. 

We proceed in a similar way to the case of even n, so we also discriminate the cases a = 0 and b = 0. In this case, we 
denote by I the inequality corresponding to (92). 

- Case a = 0. This is equivalent to b = (n + So)/2. We have that 

In,So (^ n+ 2 S0 + fi(v,So,n)d + f 0 (v,S 0 ,n), (A44) 

as in (A36). So, its second derivative is negative and we must find the minimum at the boundary; i.e., either at c = 0 
or d = 0. If we choose d = 0 (See [52] for v = 0) then we obtain I > 0 for all n, so we shall omit its discussion 
here. If we consider d = (n — So)/2, then 

In,So . n ^ 5 °) = ( 2 ) + !)(«5o + 2^ - 1), (A45) 


which is minimal for So = —2v. However, now n is odd, which implies that So must be odd as well. Hence, the 
closest odd integer values to —2v are 1 — 2i/ and —1 — 2^ which give the following coordinates (a, b, c, d): 


nil n -F 1 

0, —x- v, 0,- 


2 


2 


(A46) 
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- Case b = 0. 

Here we consider 


I n ,So (0, d) = -M 2 + g 1 (v,S Q ,n)d + g 0 (v,S 0 ,n), (A47) 

and we inspect the boundary of the domain of d. 

If d = 0, we obtain 

in,So (0,0) = Q (So + 2v + 1)(«S 0 + 2i; + 3), (A48) 


an expression which is minimal for Sq = —2(1 + v), which is an even number. As it has to be odd, the candidates 
to be considered are —3 — 2v and —1 — 2v, both of which produce a value of / = 0. The corresponding coordinates 
are 


n— 1 n n + 1 „ . 

~2 - *6 °> 0 ) . 


n — 3 n + 3 

— - v -> 0 ) 2 ^ ^ 0 ) . 


Finally, if d = (n — So)/2, we have 

r I n n - s o 

In,S 0 I U 


— (ti — 1)(S„ + 2v + 1)(25q + n(So + 2 v + l))/2, 


(A49) 


(A50) 


which is minimal for Sq = — (l + 2z/)^x§, a value which is not integer. The closest odd integers in this case are 
— (1 + 2v) and 1 — 2v, and the smallest value of I is given by the first one, for which we obtain 1 = 0 and the point 


n~ 1 n + 1 

—^ t 7 , 0 , 0 , —— 


(A51) 


Since the points (A46, A49 and A51) comprise all the vertices of P ® 11 in which J is minimum and it is 0, we deduce that 
(91) is a valid Bell inequality, as we have checked all the cases. 

□ 


Appendix B: Blocks 


In this section we discuss what is the form of Sk and Ski in the block-decomposition given by the basis (62), where k, l £ 
{x,y,z} and 

n —1 n —1 

S k = a k ] ' Skl = a k ] ® a l\ ( B1 ) 

i —0 i,j = 0 

where <j x , a y and cr z are the Pauli matrices. 

Theorem 12. Let us define m = n — 2 J and let \ £k,j) = \D% j) (g) | as in (62). Then, the J-th block pj of Sk or Ski, 

where k, l £ {x, y, z }, in the form (61 ) is given by 


{€k,j\Sy\£l,j} 

(£k,j\S z \(il,J) 
(£k,j\S xx \£,i,j) 

(£k,j\S X y\£l,j) 

(tik,j\S xz \£,i,j) 

(€k,j\Syy\€l,J ) 
{£k,j\Syz\€l,j) 

( £k,j\S zz \t;i,j) 


y/(l + 1)(2 J — 1)6 (k -l- 1) + y/(k + l)(2J-k)5(l -k- 1) 
i • sgn(/c - i)(€k,j\S x \€i,j) 

(2 J - 2k)5(k - l) 

[2k(2 J - k) - m]6(k -l) + y/(l + 1 )(l + 2)(2 J - l)(2J - l - 1)6 (k -l- 2) 

+ y/ (k + l)(fc + 2) (2 J - k)(2J -k- 1)6(1 — k — 2) 
i • sgn(fc - 1)^(1 + 1)(( + 2)(2J - l)(2J - l - 1 )S(k -1-2) 

+i • sgn (k - l)y/(k + 1 )(k + 2)(2 J - k)(2J - k - 1)6(1 - k - 2) (B2) 

(2J-1- 2l)sJ(2J -l)(l + l)S(k -l-l) + (2J-l- 2k)y/(2J - k)(k + l)S(l - k - 1) 

[2k(2J - k) - m]6(k - l)- y/(l + 1 )(l + 2)(2J - l)(2J - l - 1)6(k -1-2) 

-s/(k + l)(k + 2)(2J- k)(2J - k- 1)6(1 -k-2) 
i • sgn(fc — l)(2J — 1 — 21) \J(2J — l)(l + l)S(k — l — 1) 

+i • sgn(fc — l)(2J — 1 — 2fc)i/(2J — k)(k + 1)6(1 — k — 1) 

[(2 J — 2k) 2 — 2 J — m]S(k — l) 
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where 0 < k, l < m, 5 is the Kronecker delta function and sgn is the sign function. 


5(a) 


1 if a = 0 
0 else. 


(B3) 


Proof. We shall start by obtaining the expression of (Bl) in the computational basis. For this purpose, it will be convenient to 
introduce a bit of notation to keep the calculations simple. We will denote the elements of a 2" x 2 ra matrix A as (A)! (i-th row, 
j-th column, 0 < i,j < 2 n ), i.e., (i |x4|j). We shall also use the binary representation of i = *o*i • ■ • *n-i with i k £ {0,1} and 
the same for j, and we shall refer to i and j as words. We also define the weight of i, denoted |i|, as the number of l’s in the 
binary representation of i: |i| := \{i k : i k = 1 } |. We shall also use the bit-wise XOR function (equivalently, bitwise addition 
modulo 2) i © j := (zq © jo)(*i © j t) ■ ■ ■ (i n - 1 © jn-i)- We will say that two words i and j are at Hamming distance k if they 
differ exactly in k bits; equivalently, if |i © j| = k. 

• (Sx)j = 5(|i©j|-1). 

This expression stems from the fact that (oE)l has nonzero entries (actually, their value is 1) on those elements whose 
indices fulfill i„ © j 0 = 5 (a v) for all 0 < a < n. In other words, i and j only differ in the bit corresponding to the 
party in which a x is applied. This bit difference is a consequence of the form of a x = |0)(1| + |1)(0|. Summing over all 
parties v leads to the result. 

• ('S'y)j =i • sgn(i — © j| — 1 ). 

In this case we apply the same reasoning as in S x , but because of the form of a y = — i|0) (11 + i|l)(0|, the elements in the 
upper half are multiplied by — i whereas the elements in the lower half are multiplied by i. 

• (S z ) • = (n — 2 |i| )<5(i — j). 

This easily follows from the fact that a z is diagonal and, when we sum a z over all the parties v, those indices l with a 0 
in the u-th bit (i.e., i v = 0 ) get a +1 contribution, whereas the bits with i v = 1 get a —1 contribution. 

. (S xx )i = 2«5(|iffij|-2). 

Observe first that (cr£ l ' ) )j = <5((i © j) — (2 V + 2 W )). This is precisely the condition that bits v and w have to be 
different. Summing over all (ordered) pairs of parties v and w, the result follows. 

• (£*„)] = i • sgn(i — j)[l — (-l) 1 ^ ki ]<5(|i©j| - 2 ). 

For this case, we shall compute first the value of (cr ^&y W ' > )y Observe that (cr^ '*)j = <5((i © j) — 2 V ). and (cry^ ' 1 )j = 
—isgn(i - j)<5((iffi j) - 2 W ). Then, 

(4 V H W) )j = E(^ ) )L(4 W) ) J - = ^E^ 1 ® k ) - 2 X( k ® J) - 2”’)sgn(k-j). (B4) 

k k 

Consider the factor <5((k©j) — 2 1 ")sgn(k—j), since k and j only differ in the ic-th bit, one has that the sign of k—j is positive 
whenever j w = 0 (k w = 1) and negative whenever j w = 1 (k w = 0). Then, one can rewrite 5( (k © j) — 2 UJ )sgn(k — j) = 
^((k © j) — 2 w )[5(j w ) — 5(j w — 1)] = 5((k © j) — 2 w )(—iy w . 

(4 v) 4 W) Tj = -i(-l)'-£«(( i ® k ) - 2 t ')5((k©j) - 2 W ) = -i(-l y-(a^a^)\. (B5) 

k 

Hence, 

(4” ) 4’" ) )j = -i(-l) j ^((i©J) - (2" + 2“’)). (B 6 ) 

Now we sum (B 6 ) over all ordered pairs of particles (v, w). A convenient way to perform this sum is to join term (v. w ) 
with (w, v): 


(S*v)j = E (^^)j = -i E + (-!) J E <J((i ® j) - + n)- CB7) 

v^w 0 <v<w<n 

The term (—1) J “’ + (—l)- 3 ” can take 3 different values. It will be 0 whenever j v j w , it will be 2 whenever j v = j w = 0 
and it will be —2 whenever j v = j w = 1. Observe that only terms for which i v f j v and i w f j w are counted in the sum 
(B7). Hence, one can substitute (B7) by the following expression: 


(S X y )] = -i E f 1 - (- 1 ) 


m—m 


0 <~v<Lw<n 


sgn(i-j)5((i©j)-(2“ + 2“')) = -- [l - (-1) 


i|-|j 


sgn(i-j)(5 xx )-. (B8) 


Following the same argument as for (S xx )^, the result follows. 
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(Sxz) j = (n- |i| - |j|)<5(|i©j| - 1). 

In this case, we use the fact that S xz = S x S z — YY=o {cr x cr z )( v \ Note that 

' 0 -1 \ 


O'xO'z — 


1 0 


(B9) 


which means that XX=o ( a x a z)^ behaves essentially like S x but with a —1 instead of a 1 on the upper diagonals. Hence, 


we can write 


f n— 1 


^{vxVzY = sgn(i-j)5(|i©j| - 1) = (|i| - |j|D(|i©j| - !)■ 


(BIO) 




The last equality stems from the fact that, given a j, a number i greater than j, but at Hamming distance 1 from this j can 
only be constructed by changing one of the 0 bits of j to 1, hence ji| = |j| + 1. The opposite happens when we want to get 
a lower number: we turn one of the 1 bits of j to 0, hence lowering its weight by 1. At this point we just need to calculate 


(S xz )\ = [(^)L(^)j c ] - (|i| - j| -!) 

= £ [<5(|i © k| - l)(n- 2|k|)«y(k - j)] - (|i| - |j|)<S(|i© j| - 1) 

= (n- 2|j|)«S(|i© j| - 1) - (|i| - |j|)<S(|i© j| - 1) = (n- |i| - |j|)S(|i© j| - 1). 

(S w )j = 2(-l)^(|i©j|-2). 

In this case, we proceed in a similar way as in (S xy )y 

O v )f7 v W) )\ = ^(SyYkiSy)* = - © k - 2 u )(5(k © j - 2“)sgn(i - k)sgn(k - j) 

k k 

= ^ 5(i © k — 2 l ’)<5(k © j — 2 U ')(—1)* ,,+J “ = (-1 


(Bll) 


(B12) 


Summing (B12) over all parties, we obtain the result, since (—1 i s equal to 1 whenever i v = j w and to —1 whenever 

i v ^ j w . The Kronecker deltas in (B12) enable us to substitute (—1 by (—1) 2 Jl in this case. 


• (S yz ) ] =i-sgn(i-j)(5 xz )l. 


In this case, it suffices to change the signs the upper half of the matrix and multiply by i, as we did to go to (S v )j from 

(S zz ) ■ = [n{n — 1) — 4|i|(n — |i|)]5(i-j). 

For this last case, we apply that S zz = S Z S Z — J2v=o z&z)^■ 


(5»*)j = £( S *)k( 5 *)j k - 




= E [( n - 2 l k l) 2 ^(i - k )^( k - J)] - - j)- (B13) 


Because <5(i — k)<5(k — j) = <5(i — j)d(k — j), we obtain the result: 

{Szz)] =<5(i-j) [(n — 2|i|) 2 — n] . 


(B14) 


Let us now express | D^j)\xjj )® m / 2 in the computational basis. It will be convenient to denote A(i') = YYoYq 1 ^(* 2 a ~ 
® 2 a+i ~ 1)- Let us also denote by H = 10101010...10 and by i' • H the scalar product 1 ® 2 a+i- Then, we have that: 

2 2 - 7 - 1 /o r\- 1/2 2 m— 1 

\°2\j) = E J <Kiii - m, 1 rr m/2 = E 2- m/4 (-i) i ■ K A(i / )ii / ). 

i —0 ' ' i '—0 


(B15) 
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The matrix for the change of basis, denoted p, is then 

2 2J -12 m -l 2,7 . 

»= E E E-fe== 

i—0 i'=0 j=0 A /( 2 /)2 m / 2 

Observe thatp T p = l 2 j+i, whereas pp T is the projector onto the space Hj in (59). We shall denote |i) = |i)|i 7 ) and i = 2 m i+i 7 
in order to make the expressions more compact. 


S(\i\- j)A(i')\i)\i')(D k 2J \(^-\^. fB 16) 


(p t s xP )1 


(p T SxP)t = E^^M 


u 


= £ Ej==^a ( i')aa')«(|!® ]| - i)«(|i| - fc)J(|j| - ()- 

Using the fact that |iffi j| = |i © j| + |i' ffi j 7 |, it is convenient to distinguish the two possible cases for which <5(|i© j| — 1) 
is nonzero: |i © j| =0 and |i © jj = 1. Equivalently, i = j and i' = j 7 . We can write <5(|i © j| — 1) = 5(i — j)<5(|i 7 ffi j'| — 
1) + <5(|i ffi j| — 1)f)( i 7 — j 7 ). The first term does not lead to any contribution, because if i 7 and j 7 differ exactly in one bit, 
then either |i'| or |j'| cannot be m/2, and A(i')A(j 7 ) neutralizes any index with weight other than m/2. Hence, we are 
left with the last summand, which is 

(P T SxP)i = Y^P\(S X ) jpf 

iJ 

= £ M)1 '^ 1 ^W )A(mii®ji - w -j'Mdii - m\j\ 0 


= E 


•4 v^EYHY) 


•j.*' vWYH T 


= E 


i,j 


= E 


^RYHY) 


: A(i 7 )A(|i ffi j| - l)5(|i| - k)6{\j\ - l) 
<5(|iffij| - l)5(|i| - fc)5(|j| - Z)E A ( i7 ) 


u \/(Y) (Y) 


<5(|iffij| - l)5(|i| - fc)<5(|j| - l). 


Now we use the following fact: 

2 2J — 1 


E <KIJI - l W e j| - 1) = |i|*(|i| - l - 1) + ( 2 J - \m\i\ -1 + 1), 

j=o 


(B17) 


which is explained as follows: The left hand side of (B17) counts, among all words j, which of them have weight l and are 
at Hamming distance from a given i (they differ exactly in one bit). There are |i| of them which have weight l = |i| — 1 
(i.e., we turn a 1 into a 0) and there are 2 J — |ij of them which have weight l = |i| + 1 (i.e., we turn a 0 into a 1). 

(p T S x p)f = y Z 1 ' 1 ~ (fcj(fc - i - 1 ) + (2J - k)S(k -1+1)) 

> /©KY) 

_ k(^)S(k-l-l) + (2J-k)(^)S(k-l+l) 

men 

= (i + 1) -l-l) + (2J- QyJffi/igJSd -k- 1) 

= y/(l + l)(2J-l)S(k -l- 1) + y/(k + l)(2J-k)S{l - jfe - 1). 


• For ( p T S y p)i , the same reasoning as (p T S x p)f follows, with the additional fact that the term sgn(i — j) turns into 
sgn(k — l) because 

sgn(i-j)5(i' - j')<K|iffi j| - l)«5(|i| - k)8( |j| -1) = sgn(i - j)5(|iffi j| - l)<5(|i| - k)S(\j\ - l) = k - l, 
and (jfe — l)[S(k — l — 1) + 6(1 — k — 1)] = sgn (k — l ). 
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Let us calculate (p T S z p)^. 


( P T SzV)i = XMo 5 *)M 

IJ 


= X! (U mi A(i / )A(j / )^(i - j)5(|i| - A:)5(|j| - Z) 


IJ 


V 2m ( 2 /)( 2J ' 


= 13 ^-^===— A (i , )^(|i| - fc )^(|i| - 0 


V 2m ( 2 /)( 


2 (2 J-j 


\ - <^(|i| — Zc)i5(fc - Z) ( n ~ 2 |i| ~ 2|i'|) A 

' / /O 7\ /O 7\ /r '™ 


(VHY) “ ^ 

(Y)^ k — l)\^(n-2k- 2m/2) 


E 


' 2 J\ / 2 J\ ^ a/ 2 ™ 


(Y )( 2 

= w-oE A(i ' } = (2 J ■ 2fc)<5(fc -*)• 


We have used the fact that any function / of [i'|, when multiplied by A(i'), is a function of m/2, because only words with 
weight m/2 survive A(i'). 

Hence, we obtain ( p T S z p )f = 2(J — k)6(k — Z). 


• Let us now move on to (p T S xx p)/. 

( v T s xx p )f = 53 4 (' S '**)M 

iJ 

= 2 E ^=^L^A(i')A(j')J(|l® j| - 2)«(|i| - *)J(U| - 1 ). 

if V 2 m (Y)(Y) 


Similarly to the case of <S X , now we use the fact that <5(|i © j| — 2) = <5(i — j)<5(|i' ® j'| — 2) + <$(|i ® j| — 1 )<5(|i' ® j'| — 
1) + <5(|i® j| — 2)<5(i / — j'). The term in the middle is irrelevant, as the delta factors A(i')A(j') will cancel it out (because 
of the term <5(|i' ® j'| — 1), not both of them can have weight m/2). Let us consider the remaining terms: 


- For the case that i = j, its contribution is 


2 E 4==^A(i')A(j'W|i' © j'l - 2)<S(|i| - k)8{\i\ - l) 

w ^ m (Y)(Y) 

= 2S ( k o 53 ( ~ iy ^ 1)J - 2) 

i'.j' V 

= 2 s(k-o 53 ^^(|i , L ®j , L i-i) = -2<5(fc-z) 53 E^<y(|i , L ©j , L i-i) 

i LjL=° i LJL=° 

2 m/2 -l /9 

= -2 S(k-l) 53 = -mS{k-l), 

ij =0 V 


where we have used the encoding for the logical qubits |0 l) = 101) and jl^) = 110). Observe that A(i'')A(j / )<5(|i / ’© 
j'| — 2) = (5(|® Jl| — 1), because a pair 01 has to go to a pair 10 in order to be considered. Observe as well that 
Si' <5 (IIl ® JlI — 1 ) = m/2, as there are m/2 logical bits to flip. 

Jl 
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- And for the case i' = j', its contribution is 

2 E ( - | 1) "* ,( - 1) ' '* , A (i'){(|ia;j| - 2 ){(|i| - *)S(U| -1) 


5 ? v' 2 ”©)©' 


= 2 E 




.(5(|i©j| - 2)<5(|i| - fc)<5(|j| - l) 


\fWWD*i 


£*( 1 * 1 -*) 


5(|i|-2-0 + |i|(2J-|i|)<y(|i|-0 


2 J - |i| 

2 


A(|i| + 2 - Z) 


( k )5{k - 2 - 0 + k(2J - k)S(k -l) + { 2J ~ k )S(k + 2-1) 


2 ( 2 fc J ) 

" 7(X) 

= 2 fc (2 J k)S(k -l) + 2 ( 2 ) yJ( 2 £)/( 2 ?)S(k -2-1) + 2 ( 2J ~ k ) ( 2 k J ) / ( 2 ?)S(k + 2-1) 
= 2k(2J - k)S(k -l) + y/(l + 2 ){l + 1)(2 J - i)(2J - l - 1)6 (k -2-1) 

+ ^(2 J - k){2J - k - 1 ){k + 1 ){k + 2)6(1 -k- 2). 


We have used that, given a word i, the words j at Hamming distance 2 from i can be classified into 3 subsets: (©) 
of them having weight |j| = |i| — 2 ; (^)( 2J 77 °f them having weight |j| = |i|; and ( 2J “of them having 
weight |j| = |i| + 2. One can check that the relation ( 2 2 J ) = (^) + (^) ( 2J ^^) + ( 2J “^) holds, as it counts all the 
possibilities (all the ways to pick the 2 different bits between i and j). 


• The case (p T S xy p)^ is similar to (p T S xx p )^, so we shall point out only their differences: 

{p T SxyP) k = £?4( 5 W)jj4 

1 ’ A(i / )A(j / )sgn(i — 

2 m (Y)(Y) 

Again, we divide into the cases where i = j or i' = j'. 



j) (l - (- 1 )^) 5(|i©j| ^ 2 )J(|i| - k)6{\j\- l). 


- Case i = j and ji' © j'| =2: 

The sign function is then determined by the part corresponding to the singlet: sgn(i' — j'), and we have 




(—i) i ' N (—i)j' N 


w V 2 m (Y)(Y) 


-A(i')A(j')sgn(i'-j') l-(-l) 


i'l-IJ'l 


i(|i'©j'|- 2 )f(|i|-fc)i(|i|-l). 


But if i' and j' differ exactly in two bits, because of the factors A(i') A(j'), we only have to consider when a pair 01 
turns into a pair 10 and viceversa. So, |i'| = |j' and the whole expression is zero. 

- Case i' = j' and |i © j| = 2: Now the sign function is determined by the part corresponding to the Dicke states: 
sgn(i — j), so we have 


\2i'-K 


i£ ^===fr= A ( i ')sg n (i-j) [l - (-I)' 1 ' 2 ''' <5(|i©j| - 2)<5(|i| - k)S( |j| - l) 


= >£ 


& ]/ 2m ( 2 k J )(T 

\—> 1 


i -j \JCk )( 


= >£ 


2 + (U\ 
1 

1 

* v(YKY) 


sgn(i-j) l-(-l)' 2 J 5(|i©j| — 2)5(|i| - k)6(\j\ - 1)^2 


A(iQ 

2 m /2 


■sgn(i-j) l-(-l)^ <5(|i © j | — 2)<5(|i| — fc)5(|j| — l). 


We observe now that there are only 3 possibilities: |i| = |j| ± 2 and |i| = |j|, but this last one clearly leads to 0. 
If |i| = |j| + 2 , then i > j, because we must change exactly two bits, and they must be two ones turning into 
two zeroes, as we have the < 5 (|i © j | — 2) factor. Similarly, if |i| = |j| — 2, then i < j. We can simplify then 
sgn(i — j) = (|i| — |j|)/2. Hence, we can write the following: 


2 i E -1 ==^® JI - Wl - *)^(UI - <)• 
2 
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Now, by following the same discussion as in (p T S xx p)i 9 we arrive at the expression 


i 


k—i 

2 


2k(2J - k)S{k - l) + v/(Z + 2)(Z + 1)(2 J - Z)(2 J - l - 1)6 (k -1-2) 
+ y/( 2 ^ - fc)( 2 ^ - fc - l)(fc + 2)(fc + 1)<5(7 - k - 2)1 . 


The first term vanishes and the factor (k — l)/2 can be substituted by sgnffc — l) because of the Kronecker deltas, 
leading to the result. 

• Let us now turn to (p T S xz p)^. 


(p T S xz p)f 




E 




V 2m (Y)( 


2 J\ (2J^ 


-A(i')A(j ')(n 


HI - IllMdi© j| - i)tf(|i| - fc)«y(|j| - Z). 


As in (p T Sxp)^, it is convenient to distinguish <5(|i ® j| — 1) into the cases |i| = |j| ± 1, i' = j' or i = j, |i'| = |j'| ± 1. 
The second one is neutralized by the functions A(i') A(j'), so we only need to consider c)(|i ® j| — l)5(i / — j'): Summing 
over i' we obtain 


^ ,— 1 — —[n - |i| - |j| - 2 (m/ 2 )]^(|iej| - I j.)l i| - /.:KK|.i| - /) 

Wf) 

Recall that, given i, there are |i| words j having weight |i| — 1 and 2 J — |ij words j having weight |i| + 1. Thus, summing 
over j we obtain 

C P T S xz p )? = E 7—T-- ( 2J - k Wl - k) [|i|J(|i| - l - 1) + (2J - |i|)*(|i| -1 + 1)] 

• a/(Y)(Y) 

= (Y) ^ ( y )(y) ( 2 J - * - 0 l k6 (k - 1 - 1) + (2 J - k)5(k - l + 1)] 

= (2 J - 21 - l)y/(2J-l)(l + l)5(k - l - 1) + (2J - 2k - 1) i/(2 J — k)(k + 1)6(1 -k-1). 


Now we analyze [p T S yy p)^. 

(p T SyyP)l = E^^M 


= 2 X ^L^== a (i') a ff)(-l) 1 ^ { (l!®il - 2 M(I>I - will - «■ 


2”(V)(Y' 


This case again looks very similar to ( p T S xx p)f. Recall that £(|i® j| — 2) = A(i — j)<5(|i' ® j'| — 2) + <5(|i® j| — l)<5(|i' ® 
j'| — l) + <$(|i®j| — 2)<5(i' — j') and that the term in the middle is irrelevant, as the delta factors A(i')A(j') cancel it out. 
Hence, we distinguish the two cases: 


- Case i = j: For this case, its contribution is 


2 E ( 1)1 ^ 1)J,K A(iOA(jO(-l)^^(|i / ® j'| - 2)<5(|i| - k)5( |i| - l) 

fa V 2m (Y)(E . ; ^ 

= 2 5(k - 0 E ^ 1 4^A(iOA(j')(-l)^5(|i' ® j'| - 2). 

i'J' V 


Observe that (—1)(* + j k* = (—1)(* ®-> and that (—1)(* ®j A(i / )A(j / )<5(|i / ® j'| — 2) is nonzero only in those 
terms for which i' ® j' = 0000 • • • 001100 • • • 0000, which means that (i 7 ® j') • K = 1. Also, the combination 



43 


A(i')A(jW©j'|-2) 


only allows for words with equal weight to be considered. Hence, 

2 S(k - i) E ^A(i')A(j')(-l)^<5(|i' © j'| - 2) 

i'J' V 

= -2 6(k - 0 E i A(i')A(j')^(|i' © j'l - 2) 

i'J' V 

= -25(*-oE 

'l'Jl 

2 m/2 -l 

= -25{k-l) E ;^p=-m5(fc-0- 


- Case i' = j': For this case, it will be convenient to consider 

E(—1)^<^(|i©j| - 2)<5 (|j| - l). 

j 

There are (^) words j at Hamming distance 2 from a given i and with weight |j| = |i| — 2, which are counted 
with the factor (— l)(l 1 Hj|)/ 2 = — 1 ; there are |i|(2J — |i|) words with weight |j| = |i| which are counted with 
a factor (— l)(l 1 Hj|)/ 2 = i ; and there are ( 2J ~^) words with weight |j| = i + 2, hence counted with a factor 

(_l)(|i|-|j|)/2 = 

Now, the same argument we used for (p T S xx p)f applies, and we get a contribution from this term equal to 


2fc(2 J - k)6(k -l)- ^{l + 2 ){l + 1)(2J - 1){2J - l - 1)5 (k -2-1) 

-sJ{2J - k)(2J - k - 1 ){k + 1 ){k + 2)5(1 -k- 2). 

Summing the two contributions, we arrive at the result. 

• For the case (p T S yz p)f, we use the same argument than that for (p T S xz p)\ and, similarly to the (p T S y p)f case, we gain 
the extra factor isgn(fc — l). 

• Finally, for the case (p T S zz p)f, we have to consider 


( P T SzzP)i 


= E^( S *4 


J 

Pi 


1 J 


= E ( 1) " l<( 1>J " A(i')A(j')[(n - 2 |i |) 2 - »]4(i - j)4(|i| - fe)«(|J| - i) 

ir 

= E 1 ~, 1) " l< . < : 1 .. ) ‘ A (i')A(l')[(n - 2|i|) 2 - n]4(|i| - t)4(|i| - i) 




, 1 [(n - 2|i| - 2m/2) 2 - n]4(|i| -i)V^2 

sfClWI) 


= s(k-i)Y 

= *(* - 0 £ /— l — — i( 2J - 2 i i » 2 - 

‘/(X) 

= 6(k - l)[(2J - 2k) 2 -n] = [(2 J - 2k) 2 -2 J- m]S(k - l). 


□ 


Appendix C: Experimental errors and 2-body Bell inequalities 

Here we discuss the role of experimental errors and imperfections in detection of nonlocality with our two-body Bell in¬ 
equalities by measuring the collective observables. Our analysis is motivated by the recent experiment aiming at detection of 
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entanglement in the Dicke states [56]. We analyse how errors in the quantities Sg affect the ability to observe quantum violation 
with our Bell inequalities. 

In Sec. V E we have shown how our permutationally invariant Bell inequalities in the case when all parties perform the same 
pair of measurements can be tested with the aid of total spin components S z (i = x, y, z) and their combinations m • S along m. 
Here we take a bit more simplified approach and express our Bell inequalities in terms of only S z , S x , S%, S) and S^u, where 

s n/ 4 = (S x + S z )/v 2. 

To be more precise, let us recall that any permutationally invariant two-body Bell inequality can be written as 

Pc + o^o + PS i + + &Soi + -<Sn — (Cl) 

where Sk and Ski are defined in Eqs. (22) and (23), and at each site we use the same pair of measurements of the form 
= cos(0fc)cri^ + sin(0k)(Tx > . For simplicity, in what follows, we denote 9 0 = ip,9\ = 9. Denoting then by S u and S uv 
with u, v = x. z the permutationally invariant expectation values Si, and Su in which the measurements are taken to be simply 
<j z and a x , we can now rewrite (Cl) as 

p c + AS Z + A'S x + — S zz + DS XZ + —S xx > 0, (C2) 

where 

A = a cosip + pcosO, A'= asioip + psuup, 

B = cos 2 ip+ 28cosipcosO + ecos 2 0, 

C = 7 sin 2 ip + 28 sin<p sind + esin 2 9, 

D = 7 cos <psin ip + 8 cos </?sin 9 + 8 cos 9 sin<p + e cos 9 sin0. 

We can clearly express the quantities S u and S uv in terms of the total spin components S, as S u = 2 (S u ) and S uu = 4,(S 2 ) — n 
with u. v = x, z, and 

5,, = 2 ({S x , S z }) = 2 (2 (S 2 n/4 ) - (Si) - (S 2 Z )) . (C3) 

This allows us to rewrite (C2) in the following form 

Pc + 2 A(S Z ) + 2 A'(S X ) + |(4 (Si) - n) + ^(4 (S 2 X ) - n) + 2 D( 2 (^ /4 ) - (S 2 X ) - ( S 2 Z )) > 0, (C4) 

which can be further rewritten as: 

Pc + 2A(S Z ) + 2A'(S X ) + 2(B - D)(S 2 Z ) + 2(C - D)(S 2 X ) + 4 D(S$ /A ) - > 0. (C5) 

Having all this we can now study how the experimental errors influence the violation of our Bell inequalities. Inspired by Ref. 
[56], we consider errors of the form 


($8 )measured K A rj(S e } , (C 6 ) 

where 0 < r/ < 1 is some visibility parameter, k is an offset that depends on n, and 9 = x. z. tt/ 4. The parameter k is strongly 
related to the detection noise (i.e., the atom counting uncertainty); thus it adds independently of 77 . In this way, (C 6 ) captures 
the behavior of the experimental data of Ref. [56], in which 77 « 0.8 and n « 100: For (Si), which is ideally 0, k dominates, 
whereas for (S^), which is 0(n 2 ), 77 dominates. For (S^/ 4 ) the results of this measurement are not available, but we expect that 
both effects are taken into account reasonably by (C 6 ). 


1. Half-filled Dicke state 

In the case of Dl',/ 2 ) with even n, the ideal values of the total spin components and its moments are 

• (S x ) = (S z ) = < 5 2 ) = 0 , 

• {Sx) = n (n + 2 )/ 8 , 
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FIG. 6. On the vertical axis, we have the unnormalized quantum violation (the value of the left-hand side of (Cl)) for n = 8000; on the 
horizontal axis there is k. The lines are for r/ £ {1, 0.95,0.9,. . .} starting at rj = 1 for the blue line. 


• ( S l/i) = n (n + 2)/16. 

Recall that the Bell inequality that detects the Dicke state if,'/ 2 ) is 


n + 2 
2 


1 fn 

2 I 2 


S 00 + -Soi — 2*^11 > 0 . 


In Fig. 6 one can see the effect that errors have on the violation of Ineq. (C7). 


(C7) 


2. State for inequality (6) 

In this other case, for the Gaussian superposition of Dicke states, we have that for large n the ideal values of the total spin 
components and its moments are 

• (S g ) = 1 / 2 ^ 3 , 

• (S x ) = n/2, 

• (S 2 Z ) = 0 , 

. (S 2 )=n 2 / 4, 

• ( S l/i) = n ( n / 2 + 1 /V / 3)/4. 

Recall that the Bell inequality that is violated by these superpositions is 

2 n — 2 Sq + 2^00 — ^oi + 2 '^ 11 — ® (C 8 ) 

In Fig. 7 one can see the effect that errors have on these quantities. This case looks much more promising, as the error in (S 2 ) it 
can support may be almost an order of magnitude larger than what was achieved in the experiment [56]. 


Few remarks are in order: 

• It seems that the considerable difference in the performance with k stems from the fact that in inequality (C 8 ) the coeffi¬ 
cients are balanced, whereas in the inequality for the Dicke states (Cl) they are non-linear functions of n. Observe that 
in both cases (S 2 }j ( j ea i = 0. 

• In inequality (C 8 ), an interesting fact occurs when we take different values of the parameter 77 for (S 2 ) and (S 2 ^). Then 
one observes that the two terms in n 2 no longer cancel each other out; hence in this way, one can mimic a huge violation 
of the order of n 2 by taking one 77 larger than the other. 
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FIG. 7. On the vertical axis, we have the unnormalized quantum violation (the value of the left-hand side of Ineq. (C8)) for n = 8000; on the 
horizontal axis there is k. The lines are for r/ £ {1, 0.95,0.9,. . .} starting at rj = 1 for the blue line. 


• It seems that in order to find a Bell inequality for ffn' 2 ) that can tolerate a greater k, it would be desirable to test other 
Bell inequalities that are also violated by I)"/' 2 ), but may have more balanced coefficients; the price to pay, of course, 
will be that their effective quantum violation will be smaller. 

• Last but not least, going to smaller but controllable numbers of atoms (say 100, or so) is, of course, another way to detect 
(maybe not so ‘many’, but still) many-body nonlocality. 




